Next Article in Journal
Fast Underwater Optical Beacon Finding and High Accuracy Visual Ranging Method Based on Deep Learning
Next Article in Special Issue
Matching-Based Resource Allocation for Satellite–Ground Network
Previous Article in Journal
ACM: Accuracy-Aware Collaborative Monitoring for Software-Defined Network-Wide Measurement
Previous Article in Special Issue
Security Performance Analysis of LEO Satellite Constellation Networks under DDoS Attack
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Load Balancing Scheme for Satellite IoT Networks Based on Spatial–Temporal Distribution of Users and Advanced Genetic Algorithms

1
The School of Electronic Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China
2
Science and Technology on Communication Networks Laboratory, Shijiazhuang 050081, China
3
School of Information and Communication Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China
4
Academy for Network & Communications of CETC, Shijiazhuang 050081, China
5
China Academy of Information and Communications Technology, Beijing 100191, China
6
China Academy of Space Technology, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(20), 7930; https://doi.org/10.3390/s22207930
Submission received: 31 August 2022 / Revised: 29 September 2022 / Accepted: 12 October 2022 / Published: 18 October 2022
(This article belongs to the Special Issue Satellite Based IoT Networks for Emerging Applications)

Abstract

:
Satellite IoT networks (S-IoT-N), which have been a hot issue regarding the next generation of communication, are quite important for the coming era of digital twins and the metaverse because of their performance in sensing and monitoring anywhere, anytime, and anyway, in more dimensions. However, this will cause communication links to face greater traffic loads. Satellite internet networks (SIN) are considered the most possible evolution road, possessing characteristics of many satellites, such as low earth orbit (LEO), the Ku/Ka frequency, and a high data rate. Existing research on load balancing schemes for satellite networks cannot solve the problems of low efficiency under conditions of extremely non-uniform distribution of users (DoU) and dynamic density variances. Therefore, this paper proposes a novel load balancing scheme of adjacent beams for S-IoT-N based on the modeling of spatial–temporal DoU and advanced GA. In our scheme, the PDF of the DoU in the direction of movement of the SSP’s trajectory was modeled first, which provided a multi-directional constraint for the non-uniform distribution of users in S-IoT-N. Fully considering the prior periodicity of satellite movement and the similarity of DoU in different areas, we proposed an adaptive inheritance iteration to optimize the crossover factor and mutation factor for GA for the first time. Based on the proposed improved GA, we obtained the optimal scheme of load balancing under the conditions of the adaptation from the local balancing scheme to global balancing, and a selection of Ser-Beams to access. Finally, the simulations show that the proposed method can improve the average throughput by 3% under specific conditions and improve processing efficiency by 30% on average.

1. Introduction

Sensing and monitoring anywhere, anytime, and anyway in more dimensions are quite important for the coming era of digital twins and the metaverse [1]. The internet of things (IoT) is a network that provides the methods and infrastructures needed to achieve states sensing of object and relative information aggregation. The essential tendency of IoT in the future is global coverage and remote controlling with a higher data rate and access amount. Therefore, the IoT on the ground combined with the satellite network would be the most potent approach to greatly enlarging the scenarios of the IoT [2]. With the recent rapid developments in satellite networks, satellite IoT networks (S-IoT-N) have been a topic of general interest in the monitoring of mining, ocean shipping, and power transmission.
The key bottleneck of the above urgent application is the enhancement of the satellite network. One of the scenarios of the sixth generation (6G) is the space–air–ground integrated network (SAGIN). The Third Generation Partnership Project (3GPP) have proposed non-terrestrial networks (NTN) projects to promote the enhancement of satellite networks since 2019 [3]. Satellite internet networks (SIN) are regarded as the most possible evolution road, possessing characteristics of many satellites, low earth orbit (LEO), the Ku/Ka frequency, and a high data rate [4].
More advanced potential technologies are introduced in SIN evolution. There are some important tendencies:
(a) Satellite beams (Sat-Beam) steering in a permanent position: SIN satellites have multi-Sat-Beams, which move at high speed relative to the earth. The footprint of Sat-Beams (F-Sat-Beam) are planned by constellation [5]. During the service of the observed satellite, its Sat-Beams steer in a permanent position.
(b) Intra-serving beam (Ser-Beam)-hopping (Ser-BH): For the exploitation of a satellite-based phased array antenna (PAA), the F-SB can be further divided into different small footprints of serving beams (F-Ser-Beam) [6]. They can serve users at different places and in different periods through Ser-BH, which achieves occasion division multiplexing.
(c) More dimensions of radio resource (RR): The introduction of combining the fifth generation (5G) new radio (NR) technologies, orthogonal frequency division multiplexing access (OFDM), and non-orthogonal multiple access (NOMA) would be exploited by SIN. Using the coverage of Sat-Beams, the users can be served by the multiplex with time, frequency, power, space, and occasion divisions [7].
Meanwhile, the requests of users are not distributed uniformly in the space domain, which causes the density of each F-Ser-Beam to be quite different. In the scope of Sat-Beams, particularly the Ser-Beams at the edge, it is hard for users to choose the long-suitable Ser-Beams for load balancing. As the granularity of RR increases, the above problem becomes more serious. The relatively shifting position between the position of Sub-Satellite Point (SSP) and the user also makes this worse.
Figure 1 shows the problems of load balancing for Ser-BH in SIN. As Figure 1 shows, the first problem is how to improve the local scheme for load balancing of adjacent Ser-Beams to adapt to global scenarios. For large degradation of signal fading, the load balancing between adjacent Ser-Beams around the edge of Sat-Beams should be paid more attention to. However, there are thousands of possible scenarios for adjacent Ser-Beams with different DoU. The optimal scheme is not only to solve the problem in the case of local load balancing, but also different local scenarios around the world. The second problem is which Ser-Beams the users should be accessing to under the DoU of spatiotemporal non-uniform. The SSPs of an observed satellite are not permanent and not in the center of F-Ser-Beams, which are different from other satellites and have different times. The trajectories of SSPs are also different in the Ser-Beams serving during (SD). That means the probability density of power (PDP) of the satellite signal at the same position on the ground would vary with time, and the corresponding capacity of each Ser-Beam would also change.
Existing research on load balancing schemes for satellite networks focuses on Sat-Beam scenarios without intra-Ser-BH, and users with unified uniform distribution in the time and space domains [8,9,10,11,12,13,14,15]. Optimization of satellite network routes is the prevailing design for achieving load balancing in this research. However, the DoU are spatio-temporal non-uniform; none of the research can solve the problems of global load balancing based on local balancing schemes and the selection of access to Ser-Beams under changeable DoU. Considering the DoU have a serious relationship with geography, the users of S-IoT-N can be divided into different grids on the surface of the earth. After combining the grids with global population distributions, we can obtain multi-grids with different DoU densities. We can operate the load balancing according to the density variances of users under different movement directions from one grid to another grid. Though the densities of grids are different around the world, all scenarios can be taken as a limited set for density variances of users. Then, the basic load balancing schemes become the original reference for global schemes and different times, which can be taken as the gene of the genetic algorithm (GA). After modeling and improving the mutation and crossover of the gene, it can effectively solve the problem of the adaptation from a local balancing scheme to global balancing and the selection of Ser-Beams to access.
Finally, this paper proposes a novel load balancing scheme of adjacent beams for S-IoT-N based on the modeling of spatial–temporal DoU and advanced GA [16,17,18,19]. The main contributions of this paper are summarized as follows:
  • In contrast to existing research on load balancing under DoU of uniform, we are the first to improve these schemes by modeling the density variances of users under different moving directions. This can solve load balancing problems under spatio–temporal non-uniform DoU;
  • Fully considering the prior periodicity of satellite movement and the similarity of DoU in different areas, we propose the adaptive inheritance iteration to optimize the crossover factor and mutation factor for GA for the first time. This can enhance the efficiency and convergence speed of GA for S-IoT-N scenarios;
  • The Ser-BH scenario is totally new to SIN. We propose a load balancing scheme based on non-uniform spatial–temporal DoU and advanced GA, which can achieve better performance of total throughput for adjacent Ser-beams.
The rest of this paper is organized as follows: In Section 2, the relative works are introduced. In Section 3, the system model and problem formulation are described. Section 4 is devoted to the load balancing scheme of adjacent beams for S-IoT-N based on the modeling of spatial–temporal DoU and advanced GA. The simulations and analysis are discussed in Section 5. The summaries are provided in Section 6.

2. Related Works

Several works on satellite load balancing have been published in recent years. Their current include load balancing based on the distribution of users, load balancing based on satellite ephemeris topology, and load balancing based on Quality of Service (QoS).
The schemes of load balancing based on DoU dictate that the users decide which satellites to access by analyzing the density of equipment distributions [20,21,22,23]. W. Liu etc. proposed a routing algorithm based on segment routing for traffic return of LEO satellite networks [20]. They dynamically divide the surface of the earth into light and heavy load zones according to the relative position between gateways and the reverse slot. They then use a pre-balanced shortest path algorithm in the light load zone and use the minimum weight path defined by congestion index as the routing rule in the heavy load zone. J. Camino etc. proposed a method for optimizing the layout of satellite beams by applying mixed-integer linear programming (MILP) [21]. They designed two different sizes of spot beams to cover the area under non-uniform DoU. In their schemes, low user density areas would be covered with large-size beams, and high user density areas would be covered with small-size beams. MILP was exploited to optimize load balancing to achieve well-distributed traffic among the different beams. Syed Maaz Shahid et al. proposed a load balancing algorithm for a multi-RAT (radio access technology) network including an NTN and a TN [22]. They first offloaded the appropriate edge UEs of an overloaded cell to underutilized neighboring cells in TN. If there were any overloaded cells after the first step, they offloaded the delay-tolerant data flows of UEs to a satellite link.
The load balancing schemes based on satellite topology dictate that the users decide which satellites to access by analyzing which planned satellites would be on the services. P. Liu etc. proposed a load balancing routing scheme by hybrid-traffic-detour [24]. They calculated and designed the shortest path and a long-distance traffic detour path to forward packets and determine which path to use for forwarding to against link congestion [25]. The scheme tried to avoid a situation in which the routing path between a node and its neighbors were directed toward the same destination node. It chose the next hop based on satellite ephemeris prediction and traffic distribution, which increased the effect of load balancing. C. Dong et al. proposed a load balancing routing algorithm based on extended link states [26]. They keep all the satellite nodes informed of the link congestion state through an active state discovery and automatic congestion state release mechanism. All satellite nodes update the route table according to the link congestion state to achieve a balanced distribution of traffic load.
Load balancing schemes based on QoS dictate that the users decide which satellites to access by analyzing the average QoS of the satellite network. H. Cao et al. proposed a load balancing algorithm under the satellite network with hopping beams [27]. They divided the beams into heavy load beam (HLB) and light load beam (LLB) groups and offloaded the user terminals that had the highest packet loss rate in the HLB to an adjacent LLB. Then they offloaded the remaining user terminals that had the highest packet loss rate in the HLB to a non-adjacent LLB with hopping beams.
After reviewing the existing works on the load balancing of satellites, most research institutes were concerned with treating satellite networks as a pure transmission network. Their schemes mainly selected a suitable route to achieve efficient traffic offloading. Usually, they focused on the scenarios of the Sat-Beams without intra-Ser-BH, and the users with unified uniform distribution in the time and space domains. However, the DoU are spatio–temporal non-uniform, and none of this research has solved the problem of global load balancing based on a local balancing scheme and the selection of access to Ser-Beams under changeable DoU. On the other hand, artificial intelligence (AI) has been widely introduced to solve random communications issues.

3. System Model and Problem Formulations

3.1. Network Architecture

As is shown in Figure 2, the network architecture of S-IoT-N can be divided into space segments and ground segments. In the space segment, there are multi-satellites S = s 11 , s 12 , s 13 , , s i j , , s m n , which are composed of a satellite constellation in orbit of height; h . s i j denotes a satellite, where i and j are the i -th satellite in j -th orbit in the constellation. Next, the number of satellites is expressed as m n . There are intra-links l s i a j a s i b j b between one satellite s i a j a and another satellite s i b j b . The satellites provide users with communication coverage on the ground, which exploits the PAA to achieve dynamic beam-hopping. The coverage areas are further planned and configured according to the satellite constellation and the surface terrain of the earth. These consist of sets of F-Sat-Beams B = B 1 , B 2 , B 3 , , B N , where B i corresponds to the coverage area. All F-Sat-Beams cover the globe; the sum of coverages is 0 N S e B i . In the system using Ser-BH, the coverage of F-Sat-Beam B i is a virtual boundary for the serving satellite on the ground. There are several Ser-Beams b i j in each F-Sat-Beam B i . F-Ser-beams correspond to the actual position of the physical beams serving users. The number of F-Ser-Beams in each F-Sat-Beam can be different. The sum of the coverage area of F-Ser-Beam 0 M S e b i j should be larger than that of the F-Sat-Beams. For example, the OneWeb constellation is a near-polar-orbit satellite constellation comprising 720 satellites in 18 circular orbital planes at an altitude of 1200 km [15,16]. The satellite can be set as S = s 11 , s 12 , s 13 , , s i j , , s 40 18 . Each satellite has about 6 adjacent satellites and 16 identical, non-steerable, highly-elliptical F-Ser-Beams. The F-Ser-Beams in each F-Sat-Beam B i can be set as b i j   ,   j = 1 16 . If the users in the coverage area are served, the situation is assumed as b i j on ; otherwise, the situation is b i j off . In the ground segment, there are several earth stations (ES) and users ( U ). The ESs are the nodes that receive or transmit information to or from the satellites, whose functions are similar to the functions of a base station (BS). The users are the IoT terminals, the number of them in the world is assumed as N . The users are distributed in spatio–temporal non-uniform by geography. We will further demonstrate the user model and distributions in the next chapter.

3.2. User Model and Distribution

The DoU are analyzed in two scopes: the global scope, and the local scope. In the global scope, the user distributions are relatively permanent and follow global population distributions. This can provide useful prior knowledge. The densities of users vary greatly between each F-Ser-Beam. In the local scope, since the coverage of F-Sat-Beams is much larger than those of cells of mobile communication networks, there are very weak mutual influences between non-adjacent Ser-Beams in S-IoT-N. Therefore, the handovers between adjacent Ser-Beams based on received signal-power are not considered in our scenarios. We are concerned with the problem of which Ser-Beams the users should access for load balancing. The solutions were further improved to adapt to different density variances around the world. Next, we proposed to model the density variances of users under different moving directions, which is shown in Figure 3.
First, we built a two-dimensional satellite to ground grid map (SGGM). In an SGGM, a plane projection S u of the surface of the earth is generated. S u can be divided into a × b grids according to the configured satellite constellation, where a is the number of grids in the latitudinal direction, and b is the number of grids in the longitudinal direction. It can be deduced as:
a = a 0 × N s 0 N s 1 × H o 1 H o 0
b = b 0 × N s 0 N s 1 × H o 1 H o 0
where a 0 and b 0 are the number of grids in the directions of latitude and longitude in the case of the referenced satellite constellation, respectively. N s 0 and H o 0 are the satellite number and the orbital altitude of the referenced satellite constellation. N s 1 and H o 1 are the satellite number and orbit height of the observed satellite constellation. The observed grid G a , b is associated with an actual latitude and longitude, which is x G a 1 , G b 1 , y G a 1 , G b 1 , x G a 2 , G b 2 , y G a 2 , G b 2 , x G a 1 , G b 2 , y G a 1 , G b 2 , and x G a 2 , G b 1 , y G a 2 , G b 1 .
Second, the probability density functions (PDF) of DoU p 0 x , y of each grid were obtained by combining the SGGM and global IoT device distribution. The value of PDF is shown in the second part of Figure 3. If the number G a , b of observed satellite grids is greater than the number G a 0 , b 0 of referenced satellite grids, that is, a 1 > a 0 , p 1 x , y   can be exploited directly. If the number G a , b of observed satellite grids is less than the number G a 0 , b 0 of referenced satellite grids, that is, a 0 > a 1 , the densities p 1 x , y   can be gained by obtaining the average of the densities p 0 x , y of the surrounding grids. Next, the new PDF of DoU p g G a , b with different specifications of the SGGM was obtained.
Third, the PDF of DoU densities p g G a , b were further normalized. We designed 6 levels of normalized densities for p g G a , b , which are Q 1 , Q 2 , Q 3 , Q 4 , Q 5 , Q 6 ; Q 1 : p g G a , b 0 : 0.02 , Q 2 : p g G a , b 0.02 : 0.06 , Q 3 : p g G a , b 0.06 : 0.1 , Q 4 : p g G a , b 0.1 : 0.2 , Q 5 : p g G a , b 0.2 : 0.5 , and Q 6 : p g G a , b 0.5 : 1 .
Finally, the PDF variances with moving SSPs were obtained. Under the observation of SSP movement directions, the differences p g x , y in PDF during a Sat-Beam serving period T d from one grid to another grid are expressed as
p g G a , b = p g G a i , b i p g G a i + v d T d , b i + i v d T d .

3.3. Problems Formulations

In Equation (3), the PDF variances of different grids changing around the SGGM can be described by p g x , y to provide a foundation for solving the problems of global load balancing adaption. Here, the problems should be further analyzed in depth. In actual S-IoT-N scenarios, users are not distributed uniformly in different grids, especially for non-centered SSPs. As shown in Figure 4, two adjacent F-Sat-Beams B 1 and B 2 were observed. B 1 can be further divided into several F-Ser-Beams b 11 , b 12 , ,   b 1 n ; B 2 can be further divided into several F-Ser-Beams b 21 , b 22 , ,   b 2 n . The Ser-Beams b e d g e B 1 , B 2   around the edge of B 1 and B 2 are of most concern. The serving periods T n of different satellites in the satellite constellation serving on the same F-sat-Beam are different, which have different start-serving-times t s and over-serving-times t o . Figure 4 shows the trajectories of the SSPs in F-Ser-Beams B 1 and B 2 at the periods of T 1 and T 2 , which are set as l 11 @ B 1 & T 1 and l 21 @ B 2 & T 1 . Next, in the period T 1 , the center-point of b 1 e b e d g e B 1 , B 2 is closer to the SSP trajectory l 21 of B 2 than SSP trajectory l 11 of B 1 . In the period T 2 , the center-point of b 1 e is closer to the SSP trajectory l 12 of B 1 than the SSP trajectory l 22 of B 2 . Considering that the distances between the center-points of Ser-Beams and satellites are similar, the SSPs in the two Ser-Beams are similar too. Therefore, the problem of load balancing lies in determining which satellite’s Ser-Beams users should access, and not the access threshold. Following this problem, it can be analyzed that the constraints of optimal load balancing relate to the density variances of DoU with the SSP movement, the number of users with different services, and the differences in DoU of adjacent F-Ser-Beams.
During a specific serving period T i , Sat-Beam B i was observed in the analysis. The maximum capacity of a satellite was assumed as T r a m a x B i , the total bandwidths as B W i , and the DoU density as p g B i   , T i   , t , t T i 1   , T i . The Ser-Beams of B i are b B i = b i 1 , b i 2 , ,   b i n , the area of each Ser-Beam is s , and the corresponding PDF of DoU is p g b i k , k = 1   n . The average arrival probability of services requests is μ i t   , t T i 1   , T i . The average service traffic is φ i . The service traffic of Ser-Beams is,
T r a b i k = S p g b i k × μ i t × φ i   .
The adjacent Sat-Beam of B i is B j . The adjacent Ser-Beams are b m b e d g e B i , B j . The judgment coefficient function λ b m , is defined as
λ b m , B i = 1 ,   u s e r s   i n   b e a m   b m   a c c e s s   t o   s a t e l l i t e   B i 0 ,   u s e r s   i n   b e a m   b m   n o t   a c c e s s   t o   s a t e l l i t e   B i
Next, the traffic of the Sat-Beam is,
T r a B i = b i k b B i T r a b i k + b m b e d g e B i , B j λ b m ,   B i × T r a b m b i k b B i b e d g e B i , B j T r a b i k
We can also obtain the traffic of Sat-Beam B j with T r a B j .
Finally, the problem of load balancing for adjacent Ser-Beams of S-IoT-N can be formulated as:
  • Objectives:
Objective 1: achieve the maximum throughput of adjacent Sat-Beams:
MAX   T r a B i + T r a B j  
Objective 2: achieve the minimum waiting time for users of adjacent Sat-Beams:
MIN a v e r a g e T r a b i k / B W i
  • Conditions:
Condition 1: The PDF variances of DoU in the direction of SSPs moving in a Sat-Beam B i and serving time T i are observed.
p g B i   , T i   , t = p g B i , T i , t t   , t T i 1   , T i
Conditions 2: The PDF variances of DoU in the direction of SSPs moving between two adjacent Sat-Beams in one serving time T i are observed. The DoU of two adjacent F-Sat-Beams is different.
p g B i   , B j   , T i   , t = p g B i   , T i   , t t p g B j   , T i   , t t   , t T i 1   , T i
The probability of judgment coefficient function can be varied with density variances, which are,
ln P λ b m λ B i = 1 P λ b m λ B j = 1 = ε p g B i   , B j   , T i   , t
where p g B i   , B j   , T i   , t > 0 , indicating that the density variances of B i increase more than those of B j . Thus, the probability of accessing B i decreases, achieving load balancing.
Condition 3: The density variances of DoU in the direction of SSPs moving between two adjacent Sat-Beams in two serving times T i are observed. The trajectories of SSPs are different.
  • Constraints:
Constraint 1: The total traffic of users newly accessing the adjacent Sat-Beams is less than the traffic capacity of the Sat-Beams:
T r a B i T r a max B i
T r a B j T r a max B j  

4. Load Balancing Scheme Based on the Modeling of Spatial–Temporal DoU and Advanced Genetic Algorithms

4.1. The Modeling of the Solution Using the Original GA

In Section 3, the constraints of the problems of load balancing were proposed. A local optimal solution must be achieved to solve the problems of most other cases of load balancing around the world. At the same time, the local optimal solution should be further improved to adapt to the density variances of DoU in the spatial–temporal domain. In addition, the periodic movement of satellites can provide solution optimization with more efficient prior information. Based on the above considerations, the optimal solution was modeled by GA. The overall design of the original GA is shown in Table 1.
The specific modeling is described as follows:
  • Step 1: Modelling of Fitness Function:
The fundamental purpose of load balancing for S-IoT-N is to improve the performance of the throughput and QoS. In serving period T i , the throughput of all Sat-Beams B j should be maximized, and the average waiting time of all users should be minimized. The fitness function is defined as
f i t n e s s = k 1 T r a t o t a l k 2 T _ d e l a y a v e
where T r a t o t a l and T _ d e l a y a v e are the global throughput of S-IoT-N and the average waiting time of all users.
  • Step 2: Modelling of the Solution:
As is shown in Figure 5, we assume that in the observed area, there are k Sat-Beams B i j , and the corresponding PDF of DoU is p B i j . A solution A i j to the problem of load balancing, which achieves the optimization objective, is set as a single optimal solution output by GA. A i j are encoded in the form of natural numbers, which is a matrix composed of k natural numbers, and a m denotes the m -th access scheme of the m -th Ser-Beam. Under the possible values of each solution a m , the influence factors are the density variances of the user p g B i j , T i of the serving satellite and the difference in the density of adjacent satellite users p g B i j   , B i + i   j + j   , T i within the period T i instead of traditional random influence for the schemes to which Sat-Beams have access to. This can achieve a certain optimization direction in a specific scenario.
Under the above optimization direction, we have the solutions A i j 1   , A i j 2   , A i j 3   , , A i j n , which are the initial populations of GA optimization. The offspring individuals generated by the crossover mutation in each generation and the parent individuals form a new group and N individuals A i j 1   , A i j 2   , A i j 3   ,   , A i j N are selected using the tournament method to form the next generation populations. The new populations participate in the evolution of the next generation.
  • Step 3: Modelling the Genetic Cross:
According to the individual pairing principle of matching each other, two individuals with similar fitness, A i j x and A i j y , are selected. We judge whether the crossover occurs according to the crossover probability P a c r . A floating-point number is randomly generated in the interval [ 0 , 1 ] as the judgment factor α . If α P a c r , performs the crossover operation, randomly select l consecutive elements ( l < k ) in the individual A i j x , and exchange them with the elements in the same position in A i j y to generate two offspring individuals. The probability of genetic crossover is defined as follows:
P a c r = P a c r 1 , f i t n e s s f i t n e s s a v e P a c r 0 × f i t n e s s m a x f i t n e s s f i t n e s s m a x f i t n e s s a v e   , e l s e
where f i t n e s s is the fitness of individuals with higher fitness in an individual A i j x and A i j y , f i t n e s s a v e is the average fitness of the contemporary population, f i t n e s s m a x is the fitness of the best individual in the contemporary population, and P a c r 1 and P a c r 0 are probability constants.
  • Step4: Modelling Genetic Mutation:
For every individual, judge whether there is a genetic mutation according to the mutation probability P m u t . A floating-point number is randomly generated in the interval [ 0 , 1 ] as the judgment factor β . If β P m u t , the mutation operation is performed, and the mutation point is randomly selected in the individual A i j x . Next, the elements of the corresponding point of A i j x are changed according to the selection probability in the value set to generate a child individual. The probability of genetic mutation is defined as follows:
P m u t = P m u t 0 × t m a x t t m a x ,
where t m a x is the maximum evolution generation of the genetic algorithm, t is the current evolution generation, and P m u t 0 is a probability constant.

4.2. GA Optimization by Improving the Genetic Crossover and Mutation Based on Controllable Adaptation to the Scenario

Since the load balancing of S-IoT-N is continuous, we should fully consider the prior periodicity of satellite movements and the similarity of DoU in a different area. This can provide efficient prior information and adjust the GA according to the differences in load balancing performance between the current period and the previous period. This can control the GA’s adaptation to changes in different scenario parameters to achieve better evolution to obtain better performance. Therefore, this paper proposes a controllable adaptive genetic algorithm (CAGA) based on prior information on S-IoT-N. The optimization is mainly on the improvement of genetic crossover and genetic mutation.
  • Optimization of Cross Factor:
Before the service period T i of satellite S i starts, the incoming serving satellite obtained the optimal solution A i j T i 1 , A i j T i 2 , A i j T i 3 ,   , A i j T i h of the load balancing scheme for Ser-Beams in former h periods. They grouped a set A _ p r o i j with a prior high-fitness cross object sequentially, and the PDF variances p g B i   , T i δ   , t , δ = 1   h of each optimal solution in the corresponding period were also obtained.
At the beginning of the crossover operation, individuals in the contemporary population were divided into individuals whose fitness is higher than the average fitness of the population and individuals whose fitness is lower than the average fitness of the population. Here, whether to cross was judged according to the crossover probability shown in Equation (15). Individuals with higher fitness than the population average were selected as crossover objects with similar fitness according to the principle of matching each other. The individuals whose fitness was lower than the average fitness of the population selected cross objects from the prior high-fitness cross object set A _ p r o i j , where the probability of each element of A _ p r o i j being selected is determined by the differences factor of densities variances ω T i   , T i δ , δ = 1   h . The difference factor of PDF variances can be calculated as,
1 ω T i   , T i δ = T i 1 T i p g B i   , T i   , t p g B i   , T i δ   , t 2 d t , = T i 1 T i p g B i   , T i   , t t p g B i   , T i δ   , t t 2 d t = T i 1 T i p g B i   , T i   , t p g B i   , T i δ   , t t 2 d t = T i 1 T i p g B i   , T i   , T i δ   , t t 2 d t
The probability of A i j T i δ selected can be calculated:
P A i j T i δ = ω T i   , T i δ d = 1 h ω T i   , T i d
where 1 ω T i   , T i δ is the quantification of the differences in user density mutation of Ser-Beams in two periods. When 1 ω T i   , T i δ became smaller, the differences of two scenarios in two periods became fewer, and the probability that the load balancing scheme of this period was selected increased.
  • Optimization of Mutation Factor:
Before the service period T i of satellite S i starts, the incoming serving satellite obtained p g B i   , T i   , t of F-Sat-Beams and p g B j   , T i   , t of the adjacent F-Sat-Beam which could be accessed. In the set of any feasible solution A i j = a 1   , a 2   , a 3   , a k , the solution a m to access m -th F-Sat-beams included all possible F-Ser-beams in the m -th Sat-Beam. In the operation of mutation, the mutated value of a m was no longer randomly selected from a random probability set instead of a certain probability. The factor ω B i that affects probability can be calculated as
1 ω B i   = 1 T i T i 1 × T i 1 T i p g B i   , T i   , t d t
Furthermore, the probability of the selection of scheme access to B i is:
P B i = ω B i ω B i + ω B j = 1 T i T i 1 × T i 1 T i p g B i   , T i   , t d t 1 T i T i 1 [ T i 1 T i p g B i   , T i   , t d t + T i 1 T i p g B j   , T i   , t d t ]   = 1 T i T i 1 × T i 1 T i p g B i   , T i   , t d t 1 T i T i 1 [ T i 1 T i p g B i   , T i   , t d t + T i 1 T i p g B j   , T i   , t d t ] = T i 1 T i p g B i   , T i   , t d t T i 1 T i p g B i   , T i   , t + p g B j   , T i   , t d t
In Equation (20) we can conclude that the higher the PDF of DoU of F-Sat-Beams, the lower the probability of being selected. The overall evolution direction of the mutation operation is limited to optimization in the load balancing direction. In addition, with the observed periods and F-Ser-Beams varying, the evolution direction of the mutation can automatically adapt to real scenarios following the variances of the relative differences of PDF of adjacent F-Sat-Beams.

4.3. Advanced Load Balancing Scheme Based on Optimized GA

In the first two parts of this Section, we modeled the solution using the original GA and improved genetic crossover and mutation. According to this work, we designed an advanced load balancing scheme based on an optimized GA (LB-CAGA) in this part. Before the start of each satellite service period, the service satellite inputs the relevant prior information and user information of Ser-Beams and sets parameters for the GA. Next, the GA figures out the traffic load balancing scheme for the service period. The specific solution process is as follows:
The pseudocode of the algorithm is shown in Algorithm 1:
Algorithm1: Load balance based on GA
1.
Input: T r a b i j ; B W i ; P B i ; A i j T i δ , P A i j T i δ , δ = 1   h ;
2.
Output: A i j _ b e s t , f i t n e s s _ b e s t
3.
Create a population matrix P o p u l a t i o n N , k ;
4.
for i from 1 to N
5.
P o p u l a t i o n i , : rand   ( A i j )
6.
T r a t o t a l j = 1 k T r a b i j × P o p u l a t i o n i , k
7.
T _ d e l a y a v e T r a t o t a l B W i × k
8.
f i t n e s s i k 1 × T r a t o t a l k 2 T _ d e l a y a v e
9.
end for
10.
Row vectors in matrix P o p u l a t i o n N , k sorted by fitness in descending order
11.
f i t n e s s a v e = 1 N × i = 1 N   f i t n e s s i
12.
Perform crossover operations C r o s s o v e r   P o p u l a t i o n i , :
13.
Perform mutation operations M u t a t i o n P o p u l a t i o n i , :
14.
for i from 1 to end
15.
T r a t o t a l j = 1 k T r a b i j × P o p u l a t i o n i , k
16.
T _ d e l a y a v e T r a t o t a l B W i × j = 1 k P o p u l a t i o n i , k
17.
f i t n e s s k 1 × T r a t o t a l k 2 T _ d e l a y a v e
18.
end for
19.
Row vectors in matrix P o p u l a t i o n e n d , k sorted by fitness in descending order
20.
Tournament selects N rows to form a new population matrix
21.
if reach the maximum evolution generation then
22.
A i j _ b e s t P o p u l a t i o n 1 , :
23.
f i t n e s s _ b e s t f i t n e s s P o p u l a t i o n 1 , :
24.
else
25.
Go to step 11
26.
end if
The flow chart of the algorithm is shown as Figure 6:
The load balancing scheme inputs the service traffic of the Ser-Beams, satellite total bandwidths, the optimal solution of the load balancing scheme in the former h periods, the selection probability of optimal solution, and the selection probability of a single gene. The GA is initialized first, creates a population matrix P o p u l a t i o n N , k , and generates an initial population. Next, it calculates individual fitness and the average fitness of the group. Evolution starts when initialization is complete: First, it sorts the individuals in descending order of fitness. Second, crossover and mutation operations are performed. The generated new individuals join the original population to form a new group. Third, the fitness of the new individuals and the average fitness of the new group are calculated. Finally, the population of the next generation is selected by tournament. If this evolution generation reaches the maximum evolution generation, it outputs the optimal solution and its fitness. Otherwise, it enters the next generation evolution and repeats the crossover mutation operation.
The crossover module inputs population matrix P o p u l a t i o n N , k , the fitness of all the individuals, the average fitness of the population f i t n e s s a v e , the optimal solution A i j T i δ of load balancing scheme in the former h periods, and the selection probability P A i j T i δ of the optimal solution. The crossover module judges the crossover probability and crossover mode of the individuals according to individual fitness. The individual whose fitness is greater than the average fitness of the population takes the individual behind itself as the crossover object and randomly selects the crossover starting point and crossover length. After crossover, the two sub-individuals generated are inserted at the end of the population. Individuals whose fitness is less than the average fitness of the group select the optimal solution of the load balancing scheme in the former h periods as the crossover object according to the probability and randomly select the crossover starting point and crossover length. After crossover, the two sub-individuals generated are inserted at the end of the population. Finally, the population that has completed all crossover operations is outputted.
The specific process of the crossover operation is shown in Algorithm 2:
Algorithm 2: Crossover
1.
Input:  P o p u l a t i o n N , k ; f i t n e s s a v e ; f i t n e s s i ;   A i j T i δ , P A i j T i δ , δ = 1   h ;
2.
Output: P o p u l a t i o n e n d , k ;
3.
i = 1
4.
while i Ndo
5.
if f i t n e s s i f i t n e s s a v e  then
6.
if rand P a c r 0 × f i t n e s s m a x f i t n e s s f i t n e s s m a x f i t n e s s a v e  then
7.
Randomly choose a crossover starting node
8.
Randomly choose the crossover length
9.
P o p u l a t i o n i , k perform a crossover operation with P o p u l a t i o n i + 1 , k
10.
Sub-individuals are inserted at the end of the population
11.
i = i + 2
12.
else
13.
i = i + 1
14.
end if
15.
else
16.
if rand P a c r 1  then
17.
Select prior high fitness crossover objects by selection probability
18.
Randomly choose a crossover starting node
19.
Randomly choose the crossover length
20.
P o p u l a t i o n i , k perform a crossover operation with A i j T i δ
21.
Sub-individuals are inserted at the end of the population
22.
i = i + 1
23.
else
24.
i = i + 1
25.
end if
26.
end if
27.
end while
The mutation module inputs population matrix P o p u l a t i o n N , k , the current evolution generation t , the maximum evolution generation t m a x , and the probability of single gene selection P B i . The mutation module judges whether the individual mutates according to the adaptive mutation probability. The mutated individual randomly selects the mutated gene node and selects the mutated gene value according to single-gene selection probability. After mutation, the sub-individual generated is inserted at the end of the population. The population that has completed all mutation operations is outputted.
The specific process of the mutation module is shown in Algorithm 3:
Algorithm 3: Mutation
1.
Input: P o p u l a t i o n N , k ; t ; t m a x ; P B i ;
2.
Output: P o p u l a t i o n e n d , k ;
3.
i = 1
4.
while i Ndo
5.
if rand P m u t 0 × t m a x t t m a x then
6.
Randomly select gene node for mutation
7.
Select the result of mutation by probability
8.
The sub-individual is inserted at the end of the population
9.
end if
10.
i = i + 1
11.
end while

5. Simulation and Analysis

5.1. Simulation and Experimental Design

To verify the performance of the proposed scheme, we constructed a platform to simulate the load balancing of S-IoT-N. This platform consisted of four modules: the satellite simulation module, user simulation module, coverage simulation module, and control and analysis module, as per Figure 7. The satellite simulation module was mainly used to simulate satellite nodes under different satellite configurations, which affected the size of satellite signal coverage on the ground, as well as the SSPs. The user simulation module was mainly used to simulate different DoU, which related to the population distribution and service behavior. The coverage simulation module mainly simulates F-Sat-Beams, F-Ser-Beams, the trajectory of SSPs on the ground, and the constraint boundary corresponding to the population distribution and determined constellation. The control and analysis module was mainly used to complete scenario realization, parameter configuration, and algorithm selection, as well as analysis of the learning ability of the algorithm, the throughput, and the efficiency of the traffic balancing method.
Using the simulation platform, we verified the performance of the GA’s learning efficiency and load balancing.
The parameters of GAs are configured in Table 2.
The parameters of load balancing for S-IoT-N are configured in Table 3. Load balancing simulation scenario considering 9 F-Ser-Beams at the edge of 3 F-Sat-Beams. Similar parameters can be found in the article [20].

5.2. Performance of the Improved Genetic Algorithm and Analysis

The solution efficiency and optimization effect of the three genetic algorithms were simulated and compared, and the statistics of the simulation results are as follows:
It can be seen from the statistical chart that under the condition of different population sizes in Figure 8, CAGA has higher optimization solution efficiency. The reason for this is that CAGA exploited prior information to limit the convergence direction of the algorithm, which accelerated the convergence speed of GA. The PAGA and traditional GA only choose better solutions in the selection module. However, the CAGA chose better solutions in the selection module, crossover module, and mutation module, which make the CAGA more efficient. Taking a population size of 60 as an example, the CAGA figured out the global optimal solution in the 1000th generation, which is 400 generations less than the probability adaptive genetic algorithm (PAGA), and nearly 1000 generations less than the traditional GA. Under the conditions of four different population sizes, the average generation used by the CAGA to search for the optimal solution is 400 generations less than that of the PAGA on average, and 700 generations less than that of the traditional GA.

5.3. The Performance of Load Balancing for S-IoT-N and Analysis

In the above-mentioned simulation environment of satellite traffic load balancing, LB-CAGA is compared with the most widely studied integrated weighted access scheme (IWAS). The comparison of load balancing effects of the two schemes is shown in Figure 9.
As we can see from Figure 9, the system throughput of LB-CAGA and IWAS is consistent when the user density of the adjacent F-Sat-Beam is 0.5 and 0.6. When the user density of the adjacent F-Sat-Beam is 0.7, and the user density of the target F-Sat-Beam is less than 0.74, the system throughput of LB-CAGA and IWAS is consistent. When the user density of the adjacent F-Sat-Beam is 0.7, and the user density of the target F-Sat-Beam is more than 0.74, the system throughput of LB-CAGA is on average 0.4G higher than that of IWAS. When the user density of the adjacent F-Sat-Beam is 0.8, and the user density of the target F-Sat-Beam is less than 0.62, The system throughput of LB-CAGA and IWAS is consistent. When the user density of the adjacent F-Sat-Beam is 0.8, and the user density of the target F-Sat-Beam is greater than 0.62 and less than 0.82, the system throughput of LB-CAGA is on average 0.4G higher than that of IWAS. When the user density of the adjacent F-Sat-Beam is 0.8, and the user density of the target F-Sat-Beam is more than 0.82, the system throughput of LB-CAGA and IWAS is consistent.
From the above analysis, we found that LB-CAGA does not achieve more system throughput than IWAS when the user density of F-Sat-Beam is too low or too high but achieves more system throughput when the user density of the F-Sat-Beam is high and still within the capacity of the serving satellite. The reasons are as follows:
When the user density of the F-Sat-Beam is low, the traffic load is still within the capacity of the serving satellite. Both schemes allow for the serving satellite to fully accommodate all user traffic. Consequently, the system throughput of both schemes is the same as the total traffic load. When the user density of the F-Sat-Beam is too high, the traffic load exceeds the capacity of the service satellite. Neither of the two schemes can make the serving satellite bear more traffic load. As a result, the system throughput of both schemes is the same as the capacity of the serving satellite. When the user density of F-Sat-Beam is high but still within the capacity of serving satellite, LB-CAGA can improve the total throughput of the system to a certain extent and maximize the efficiency of the satellite system. Under this condition, LB-CAGA has a better load-balancing effect than IWAS.

6. Conclusions

This paper proposes a novel load balancing scheme of adjacent beams for S-IoT-N based on the modeling of spatial–temporal distribution of users and advanced GAs. The main conclusions are as follows:
(1) We modeled the PDF of DoU in the direction of movement of SSP trajectory, which provided a multi-directional constraint for non-uniform distribution users in S-IoT-N. Compared with the existing model with permanent global distribution and pure random distributions, the proposed model can characterize the PDF variances of DoU more correctly in the scenario of a highly dynamic satellite.
(2) The crossover factor and mutation factor in GAs are proposed to optimize with prior information such as the periodicity of satellite movement and the proposed model of DoU, which can better improve the efficiency of GAs and the performance of load balancing in S-IoT-N than other existing methods.
(3) Based on the proposed improved GA, we obtained the optimal load balancing scheme under the conditions of adaptation from the local balancing scheme to global balancing and the selection of Ser-Beam access. In the scenario of extremely non-uniform DoU and dynamic density variances, advances in beam-hopping were fully realized.
The proposed load balancing scheme has many application scenarios in S-IoT-N. For example, S-IoT-N can make up for the lack of terrestrial networks caused by earthquakes and other natural disasters; this scheme can effectively relieve the sudden increase of traffic load pressure and keep the communication link unblocked. Moreover, when using sensors to monitor and transmit the circuit status information through S-IoT-N, this scheme can ensure the timeliness of the information, so as to take timely measures for various emergencies.

Author Contributions

Conceptualization, W.L., Z.D., K.W. and D.W. (Dongdong Wang); methodology, W.L., Z.D., K.W. and D.W. (Dongdong Wang); software, Z.D.; validation, Z.D., W.L., K.W. and D.W. (Dongdong Wang); formal analysis, Z.D., Y.L. (Yang Liu), Y.D. and D.W. (Da Wan); investigation, Z.D.; resources, Z.D., W.L., Y.L. (Yicheng Liao) and B.X.; data curation, Z.D.; writing—original draft preparation, Z.D.; writing—review and editing, Z.D., W.L, Y.L. (Yang Liu) and G.W.; visualization, Z.D.; supervision, W.L., K.W. and D.W. (Dongdong Wang); project administration, W.L., Z.D., K.W. and D.W. (Dongdong Wang); funding acquisition, W.L., K.W. and D.W. (Dongdong Wang). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the China National Key Research and Development Plan (2022YFB2902605 and 2020YFB1808005), Science and Technology on Communication Networks Laboratory Foundation project (HHX22641X003), the Science and Technology on Communication Networks Laboratory Foundation project (FFX22641X010), the 13th 5-year plan Civil Aerospace Technology Advance Research Project (D030301), the Hebei Province High-level Talent Funding Project (B2021003032), and the New technology research university cooperation project (SKX212010010).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Munoz, R.; Vilalta, R.; Yoshikane, N.; Casellas, R.; Martinez, R.; Tsuritani, T.; Morita, I. Integration of IoT, Transport SDN, and Edge/Cloud Computing for Dynamic Distribution of IoT Analytics and Efficient Use of Network Resources. J. Light. Technol. 2018, 36, 1420–1428. [Google Scholar] [CrossRef]
  2. Bhat, J.R.; Alqahtani, S.A. 6G Ecosystem: Current Status and Future Perspective. IEEE Access 2021, 9, 43134–43167. [Google Scholar] [CrossRef]
  3. Hosseinian, M.; Choi, J.P.; Chang, S.-H.; Lee, J. Review of 5G NTN Standards Development and Technical Challenges for Satellite Integration With the 5G Network. IEEE Aerosp. Electron. Syst. Mag. 2021, 36, 22–31. [Google Scholar] [CrossRef]
  4. Letaief, K.B.; Shi, Y.; Lu, J.; Lu, J. Edge Artificial Intelligence for 6G: Vision, Enabling Technologies, and Applications. IEEE J. Sel. Areas Commun. 2022, 40, 5–36. [Google Scholar] [CrossRef]
  5. Tang, J.; Bian, D.; Li, G.; Hu, J.; Cheng, J. Resource Allocation for LEO Beam-Hopping Satellites in a Spectrum Sharing Scenario. IEEE Access 2021, 9, 56468–56478. [Google Scholar] [CrossRef]
  6. Tusha, A.; Dogan, S.; Arslan, H. A Hybrid Downlink NOMA With OFDM and OFDM-IM for Beyond 5G Wireless Networks. IEEE Signal Process. Lett. 2020, 27, 491–495. [Google Scholar] [CrossRef]
  7. Liu, Y.; Zhang, S.; Mu, X.; Ding, Z.; Schober, R.; Al-Dhahir, N.; Hossain, E.; Shen, X. Evolution of NOMA Toward Next Generation Multiple Access (NGMA) for 6G. IEEE J. Sel. Areas Commun. 2022, 40, 1037–1071. [Google Scholar] [CrossRef]
  8. Liu, S.; Lin, J.; Xu, L.; Gao, X.; Liu, L.; Jiang, L. A Dynamic Beam Shut Off Algorithm for LEO Multibeam Satellite Constellation Network. IEEE Wirel. Commun. Lett. 2020, 9, 1730–1733. [Google Scholar] [CrossRef]
  9. del Portillo, I.; Cameron, B.G.; Crawley, E.F. A technical comparison of three low earth orbit satellite constellation systems to provide global broadband. Acta Astronaut. 2019, 159, 123–135. [Google Scholar] [CrossRef]
  10. Tall, A.; Altman, Z.; Altman, E. Self-Optimizing Load Balancing With Backhaul-Constrained Radio Access Networks. IEEE Wirel. Commun. Lett. 2015, 4, 645–648. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, H.Y.; Kim, H.; Cho, Y.H.; Lee, S.-H. Self-Organizing Spectrum Breathing and User Association for Load Balancing in Wireless Networks. IEEE Trans. Wirel. Commun. 2016, 15, 3409–3421. [Google Scholar] [CrossRef]
  12. Nishiyama, H.; Kudoh, D.; Kato, N.; Kadowaki, N. Load Balancing and QoS Provisioning Based on Congestion Prediction for GEO/LEO Hybrid Satellite Networks. Proc. IEEE 2011, 99, 1998–2007. [Google Scholar] [CrossRef]
  13. Fang, X.; Feng, W.; Wei, T.; Chen, Y.; Ge, N.; Wang, C.-X. 5G Embraces Satellites for 6G Ubiquitous IoT: Basic Models for Integrated Satellite Terrestrial Networks. IEEE Internet Things J. 2021, 8, 14399–14417. [Google Scholar] [CrossRef]
  14. Wang, C.; Liu, L.; Ma, H.; Xia, D. A Joint Optimization Scheme for Hybrid MAC Layer in LEO Satellite Supported IoT. IEEE Internet Things J. 2021, 8, 11822–11833. [Google Scholar] [CrossRef]
  15. Zhang, Z.; Jiang, C.; Guo, S.; Qian, Y.; Ren, Y. Temporal Centrality-Balanced Traffic Management for Space Satellite Networks. IEEE Trans. Veh. Technol. 2017, 67, 4427–4439. [Google Scholar] [CrossRef]
  16. El-Shorbagy, M.A.; El-Refaey, A.M. Hybridization of Grasshopper Optimization Algorithm With Genetic Algorithm for Solving System of Non-Linear Equations. IEEE Access 2020, 8, 220944–220961. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Yu, W.; Chen, X.; Jiang, J. Parallel Genetic Algorithm to Extend the Lifespan of Internet of Things in 5G Networks. IEEE Access 2020, 8, 149630–149642. [Google Scholar] [CrossRef]
  18. Zhang, Q.; Yang, S.; Liu, M.; Liu, J.; Jiang, L. A New Crossover Mechanism for Genetic Algorithms for Steiner Tree Optimization. IEEE Trans. Cybern. 2020, 52, 3147–3158. [Google Scholar] [CrossRef]
  19. Li, J.; Li, L. A Hybrid Genetic Algorithm Based on Information Entropy and Game Theory. IEEE Access 2020, 8, 36602–36611. [Google Scholar] [CrossRef]
  20. Liu, W.; Tao, Y.; Liu, L. Load-Balancing Routing Algorithm Based on Segment Routing for Traffic Return in LEO Satellite Networks. IEEE Access 2019, 7, 112044–112053. [Google Scholar] [CrossRef]
  21. Camino, J.-T.; Artigues, C.; Houssin, L.; Mourgues, S. Mixed-integer linear programming for multibeam satellite systems design: Application to the beam layout optimization. In Proceedings of the 2016 Annual IEEE Systems Conference (SysCon), Orlando, FL, USA, 18–21 April 2016; pp. 1–6. [Google Scholar]
  22. Shahid, S.M.; Seyoum, Y.T.; Won, S.H.; Kwon, S. Load Balancing for 5G Integrated Satellite-Terrestrial Networks. IEEE Access 2020, 8, 132144–132156. [Google Scholar] [CrossRef]
  23. Kyrgiazos, A.; Evans, B.; Thompson, P. Irregular beam sizes and non-uniform bandwidth allocation in HTS multibeam satellite systems. In Proceedings of the 31st AIAA International Communications Satellite Systems Conference (ICSSC), Florence, Italy, 14–17 October 2013. [Google Scholar]
  24. Liu, P.; Chen, H.; Wei, S.; Li, L.; Zhu, Z. Hybrid-Traffic-Detour based load balancing for onboard routing in LEO satellite networks. China Commun. 2018, 15, 28–41. [Google Scholar] [CrossRef]
  25. Liu, J.; Luo, R.; Huang, T.; Meng, C. A Load Balancing Routing Strategy for LEO Satellite Network. IEEE Access 2020, 8, 155136–155144. [Google Scholar] [CrossRef]
  26. Dong, C.; Xu, X.; Liu, A.; Liang, X. Load balancing routing algorithm based on extended link states in LEO constellation network. China Commun. 2022, 19, 247–260. [Google Scholar] [CrossRef]
  27. Cao, H.; Su, Y.; Zhou, Y.; Hu, J. QoS Guaranteed Load Balancing in Broadband Multi-Beam Satellite Networks. In Proceedings of the 2019 IEEE International Conference on Communications (ICC), Shanghai, China, 20–24 May 2019; pp. 1–6. [Google Scholar]
Figure 1. The problems of load balancing for Ser-BH in SIN.
Figure 1. The problems of load balancing for Ser-BH in SIN.
Sensors 22 07930 g001
Figure 2. The network framework of S-IoT-N.
Figure 2. The network framework of S-IoT-N.
Sensors 22 07930 g002
Figure 3. The user model and distribution of S-IoT-N.
Figure 3. The user model and distribution of S-IoT-N.
Sensors 22 07930 g003
Figure 4. The case of different results at different times and F-Ser-Beams.
Figure 4. The case of different results at different times and F-Ser-Beams.
Sensors 22 07930 g004
Figure 5. Inheritance or mutation of genetic control considering temporal variation.
Figure 5. Inheritance or mutation of genetic control considering temporal variation.
Sensors 22 07930 g005
Figure 6. Flow Chart of load balancing method for satellite IoT beam.
Figure 6. Flow Chart of load balancing method for satellite IoT beam.
Sensors 22 07930 g006
Figure 7. The simulation platform of load balancing for S-IoT-N.
Figure 7. The simulation platform of load balancing for S-IoT-N.
Sensors 22 07930 g007
Figure 8. The performance of improved GA.
Figure 8. The performance of improved GA.
Sensors 22 07930 g008
Figure 9. The performance of improved load balancing schemes.
Figure 9. The performance of improved load balancing schemes.
Sensors 22 07930 g009
Table 1. The solution modeling of load balancing by GA.
Table 1. The solution modeling of load balancing by GA.
Genetic Algorithm ModeModeling of Satellite User Load Balancing Based on Genetic Algorithms
Fitness FunctionThe throughput of adjacent Sat-Beam B i , B j in the satellite service period T i and average waiting time for users of adjacent Sat-Beam B i , B j .
Fitness RulesThe fitness increases with the increase of throughput and decreases with the increase of the average waiting time for users.
Single Optimal SolutionThe load balancing scheme A i j of the local k Ser-Beams b i j when the user density is p b i j and the serving satellite is S i j .
Solution EncodingNatural number encoding.
Solution Code A i j = a 1   , a 2   , a 3   , a k , a m   1 m k denotes the access scheme of the m -th F-Sat-Beams.
Factors Influencing the SolutionThe densities variances p g b i j   , T i of DoU in the direction of SSPs moving in a Ser-Beam b i j and serving time T i , differences with that of adjacent Ser-Beams p g b i j   , b i + i   j + j   , T i .
Selected Set of Solutionsn load balancing schemes A i j 1   , A i j 2   , A i j 3   , , A i j n of local k Ser-Beams b i j when the user density is p b i j and the serving satellite is S i j
A Set of Solutions Selected According to FitnessThe N schemes with the highest fitness among the existing load balancing schemes A i j 1   , A i j 2   , A i j 3   ,   , A i j N
The Process of Coding CrossoverSelect two load balancing schemes A i j x   , A i j y with similar fitness according to the crossover probability. Randomly select the points for crossover and exchange the elements of the corresponding points of the two solutions.
The Process of Coding MutationSelect a load balancing scheme A i j x based on mutation probability. The mutation points are randomly selected, and the elements of the corresponding points of A i j x are changed according to the selection probability in the value set.
Table 2. Simulation parameter configuration for GA.
Table 2. Simulation parameter configuration for GA.
CAGAPAGAGA
CrossoverTwo-point crossoverTwo-point crossoverTwo-point crossover
MutationSingle-point mutationSingle-point mutationSingle-point mutation
SelectTournamentTournamentTournament
Crossover ProbabilityBase probability 0.6Base probability 0.6fixed probability 0.6
Mutation ProbabilityBase probability 0.1Base probability 0.1fixed probability 0.1
Population Size40/60/
80/100
40/60/
80/100
40/60/
80/100
Evolution Generation0:200:20000:200:20000:200:2000
Elitist Preservationuseuseuse
Termination ConditionReach the maximum evolutionary generationReach the maximum evolutionary generationReach the maximum evolutionary generation
Table 3. Simulation parameters configuration of load balancing.
Table 3. Simulation parameters configuration of load balancing.
Number of F-Ser-Beam9
Traffic Value of F-Ser-Beam400M
User Density of Target F-Sat-Beam0.5: 0.04: 0.9
User Density of Adjacent F-Sat-Beam0.5/0.6/0.7/0.8
Satellite Total Bandwidth2G
Average Traffic Value of User Request15M
Average Arrival Rate of User Request0.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, W.; Dong, Z.; Wang, K.; Wang, D.; Deng, Y.; Liao, Y.; Liu, Y.; Wan, D.; Xu, B.; Wu, G. A Novel Load Balancing Scheme for Satellite IoT Networks Based on Spatial–Temporal Distribution of Users and Advanced Genetic Algorithms. Sensors 2022, 22, 7930. https://doi.org/10.3390/s22207930

AMA Style

Lin W, Dong Z, Wang K, Wang D, Deng Y, Liao Y, Liu Y, Wan D, Xu B, Wu G. A Novel Load Balancing Scheme for Satellite IoT Networks Based on Spatial–Temporal Distribution of Users and Advanced Genetic Algorithms. Sensors. 2022; 22(20):7930. https://doi.org/10.3390/s22207930

Chicago/Turabian Style

Lin, Wenliang, Zewen Dong, Ke Wang, Dongdong Wang, Yaohua Deng, Yicheng Liao, Yang Liu, Da Wan, Bingyu Xu, and Genan Wu. 2022. "A Novel Load Balancing Scheme for Satellite IoT Networks Based on Spatial–Temporal Distribution of Users and Advanced Genetic Algorithms" Sensors 22, no. 20: 7930. https://doi.org/10.3390/s22207930

APA Style

Lin, W., Dong, Z., Wang, K., Wang, D., Deng, Y., Liao, Y., Liu, Y., Wan, D., Xu, B., & Wu, G. (2022). A Novel Load Balancing Scheme for Satellite IoT Networks Based on Spatial–Temporal Distribution of Users and Advanced Genetic Algorithms. Sensors, 22(20), 7930. https://doi.org/10.3390/s22207930

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop