随着全球通信业务的迅速增长，卫星通信系统在各领域应用愈发广泛^[1]。未来低轨卫星网络可成为人们日常生活中的重要组成部分^[2]。低轨卫星通信系统凭借其覆盖范围广、传播时延小和终端设备发射功率低等优势，使借助LEO卫星进行数据传输越发受到广泛关注^[3]。同时，星间链路的使用既满足高质量和高数据安全性等特殊应用需求，也引发了业内对卫星星座网络路由问题的持续关注，其中负载均衡问题是路由算法研究的重要部分，也是基于星间链路的卫星网络设计中的重要问题^[4]。

相较于中高轨卫星，低轨卫星具有低时延特征，理论上能够支持时延敏感型垂直应用，如实时视频传输、远程工业控制、远程医疗等。因此，端到端时延成为低轨卫星网络负载均衡路由算法设计中需要考虑的重要因素之一^[5]。

为了满足卫星通信系统星上数据就近落地需求，提高系统总吞吐量，现有的卫星通信系统(如OneWeb等)都具有分布在全球各地的地面信关站^[6]。由于全球用户的非均匀分布及其活跃度随时间动态变化，卫星系统数据流通过地面站就近下行会引发空间段拥塞并加重馈线负担，这种拥塞模式是一种漏斗型的新型拥塞模式且极易引起雪崩式拥塞，进而引起空间段端到端通信的性能劣化^[7]。

漏斗型拥塞形态随着网络状态不断动态变化，图 1以每个地面站同时跟踪4颗卫星为例给出一种特殊的漏斗型拥塞形态。如图 1所示，由于卫星S₇、S₁₁、S₁₂、S₁₆与地面站之间的馈线传输能力有限，而各个方向通过地面站下行的非均匀数据流向地面站上方卫星簇不断汇聚，由此形成一个漏斗型拥塞区域。目前的负载均衡技术并未针对此类拥塞模式进行优化。

卫星网络多径路由策略利用卫星网络具有动态可预测的拓扑形态以及卫星节点间的天然多径的特点进行设计^[8]。相对单条路径的路由策略而言，多径路由的分流策略更加灵活。CEMR(compact explicit multi-path routing)算法首先提出卫星网络中的多径路由策略，综合考虑排队时延与传播时延计算路由表，但具体的多径实现方式并未详述^[9]。Taleb等^[10]给出一种显式负载均衡(explicit load balancing, ELB)算法，该算法监控本地拥塞状态，当发生拥塞时及时通知上游将通过本地拥塞卫星的χ%的流量通过备选路径进行转发。在此基础上Song等^[11]提出交通灯智能路由策略(traffic-light based intelligent routing strategy, TLR)，采用交通灯的“红绿灯”概念将本地队列与卫星整体队列的情况分为3级，综合考虑本地与下一跳节点的状态选择转发策略，如果最优路径及备选路径均为红灯则不适合转发，此时将数据分组存入等待区，直到任一路径恢复为非红灯状态后发出。ELB算法和TLR算法虽然能在一定程度上缓解拥塞，缩短端到端时延，但其分流策略不具有全局视野，容易陷入局部最优。Liu等^[12]进一步提出一种基于混合分流策略的负载均衡路由策略(hybrid traffic detour load balancing routing，HLBR)，该策略将长程绕行与分布式分流2种方式相结合以实现高效自适应负载均衡。但HLBR算法的分流策略较为复杂，牺牲了一定的时间复杂度及空间复杂度。

综上所述，已有的路由算法都没有针对性地解决地面站就近下行引起的漏斗型拥塞问题，而此类拥塞会严重影响低轨卫星通信系统空间段端到端通信的负载均衡及通信信息的时效性。为了满足星地传输需求，保障通信时效性并实现空间段负载均衡，本文从全网链路代价计算与多径计算2个方面对TLR算法进行改进，提出一种基于区域分流的多径搜索负载均衡路由算法(regional-traffic-detour multipath search load balancing routing algorithm, RMLBR)。

1 系统模型 1.1 星座模型

本文利用中国科学院微小卫星创新研究院某在研类铱星极轨道(walker polar)星座构建星座模型并进行路由算法设计。该星座包含Num_total=Num_orbit×Num_satperorb颗卫星及若干地面信关站，其中Num_orbit代表轨道数目，Num_satperorb代表每轨卫星数目, 全网所有卫星节点构成卫星集合S={S_i|1≤i≤Num_total}。每颗卫星都具有4条节点到节点的双工星间链路，最多可连接4颗邻近卫星，其中2条为同轨链路，2条为轨间链路。当卫星经过极区上方与反向缝处时轨间链路关闭，且卫星可以与其覆盖范围内的终端设备及地面信关站建立星地链路。此外，每颗卫星各链路的发射机中均配置参数一致的缓存队列用于缓存数据分组。

由于卫星的动态运动及地球的自转，每个卫星覆盖的范围及卫星之间的连接关系持续变化。为了便于研究，本文通过虚拟拓扑法将不断运动的实际卫星一一映射为静态的虚拟卫星，并将每颗虚拟卫星与一个固定的覆盖范围进行绑定。当实际卫星运动时，其对应的虚拟卫星也会随之变化。覆盖区域的划分与卫星星座构型有关，根据星座构型将地球表面划分为Num_total个区域。

1.2 流量模型

本文使用文献[14]中的数据及方法并结合中国科学院该在研星座试运营阶段的星座构型建立流量模型作为算法的参考输入。

首先，根据该星座Walker 72/6/3π型星座构型，利用虚拟拓扑法将地球表面划分为72个30°×30°的区域，并为每一区域计算得到一个静态设备密度指数SDI_i(static device density index)如图 2所示。此外，由于全网流量分布还具有时变性，因此计算流量比例ρ_h随时间变化情况如图 3所示。

然后，为了有效缓解星地传输带来的空间段拥塞，使用文献[7]中的方法对不同目的节点的流量进行分类：将需要通过地面信关站下行并接入地面核心网的流量称为星地流量(satellite to ground traffic，SGT)，而无需经过地面站传输的流量则称为端到端流量(satellite to satellite traffic，SST)。并根据文献[6]，对不同类型的流量分别计算对应的卫星S_i与卫星S_j间流量需求指数TDI_ij(traffic demand index，TDI)，具体公式如下

$ \mathrm{TDI}_{i j}=\left\{\begin{array}{l} \frac{\left(\mathrm{SDI}_{i} \times \mathrm{SDI}_{j}\right)^{\gamma}}{d_{i j}{ }^{\delta}}, \mathrm{SGT} ,\\ \left(\mathrm{SDI}_{i} \times \mathrm{SDI}_{j}\right)^{\gamma}, \mathrm{SST}. \end{array}\right. $

(1)

其中：d_ij为2颗卫星间距离，设置星地流量系数γ=0.5，δ=2.0；设置端到端流量系数γ=0.8。

最后，文献[12]中提出流量模型的建立还受到时间的影响，因此在卫星间流量需求指数TDI_ij基础上计算卫星间实时流量需求，计算公式如下

$ \begin{gathered} \lambda_{\rho_{\mathrm{h}}, i, j}=\frac{\mathrm{TDI}_{i, j}}{\sum\nolimits_{i=1}^{\mathrm{Num}_{\text {total }}} \sum\nolimits_{j=1}^{\mathrm{Num}_{\text {total }}} \mathrm{TDI}_{i, j}} \times \rho_{h} \times \frac{A}{3600}, \\ i \neq j, \end{gathered} $

(2)

其中：A为全网全天流量总和，流量模型中单位时间产生的数据分组服从泊松分布，由此可得卫星平均数据生成率为$ \sum_{j=1, i \neq j}^{\mathrm{Num}_{\text {total }}} \lambda_{\rho_{\mathrm{h}}, i, j}$。

2 算法描述

本文提出的RMLBR算法分为全网信息收集、链路代价计算、多径计算以及多径转发策略4个部分。算法开始时先进行全网状态信息收集建立全网信息库，再根据收集得到的信息分区域计算链路代价作为多径计算的输入，然后结合约束条件进行搜索得到一条最优路径及一条备选路径，最后在转发过程中根据当前网络状态使用“红绿灯”策略选择下一跳。

2.1 全网信息收集

有效准确地收集全网状态是保障分流策略及多径计算有效性的重要措施。在全网信息收集阶段，使用文献[11]轨道发言人策略进行网络状态收集。轨道发言人策略如图 4所示，低轨卫星网络每个轨道均设置一颗发言人卫星，轨内非发言人卫星收集本地状态并发送至发言人卫星，发言人卫星收集本轨内信息后生成轨道信息包发送至其他轨道的发言人卫星，此外，发言人卫星接收其他轨道的信息包，并转发给轨内非发言人卫星，最后建立起全网信息库。

2.2 链路代价计算

合理计算全网链路代价可以有效地描绘网络中各条链路及各区域状态，从而使算法具有全局视野并提高算法性能。为了更加细致地反映全网状态以及地面站上方空间段的潜在拥塞可能，引导算法进行分流，在全网状态信息的基础上分区域计算链路代价，避免算法陷入局部最优。

本文引用文献[7]中站域的概念，将地面站上方易拥塞的空间段区域称为站域(station area，SA)，站域上方卫星集合称为站域卫星(station area satellite，SAS)，其余卫星为非站域卫星(non station area satellite，nSAS)。对不同区域分别计算链路代价以实现多径计算。

文献[7]中以站域指数SI_i来衡量当前区域受地面站就近星地传输的影响程度进而划分站域，其中站域指数SI_i的计算使用线性模型，这样虽然能够简化运算便于仿真，但难以准确地刻画各个因素与站域指数SI_i之间的关系，使得站域的划分不够精确。

为了更准确地描绘站域形态及地面站就近星地传输对各个区域带来的潜在拥塞风险，本文对站域指数SI_i的计算方法进行改进。站域指数SI_i取值受到静态设备密度指数SDI_i、卫星覆盖区域中心与地面站之间的距离SGd_i以及地面用户活跃指数UAI_i这3种因素的影响，其中距离地面站越近的区域汇聚的星地流量越多，发生拥塞的可能性就越大，因此站域指数SI_i与卫星覆盖区域中心与地面站之间的距离SGd_i成反比；而覆盖区域的静态设备越多，用户越活跃，则该区域产生的流量也相应较大，故站域指数SI_i与静态设备密度指数SDI_i及地面用户活跃指数UAI_i成正比。综上所述对站域指数SI_i建模如下：

$ \mathrm{SI}_{i}=\frac{\mathrm{SDI}_{i}{ }^{\kappa} \times \mathrm{UAI}_{i}{ }^{\lambda}}{\mathrm{SGd}_{i}{}^{\mu}}, $

(3)

$ \mathrm{UAI}_{i}=\frac{\rho_{\mathrm{h}} \times \mathrm{SDI}_{i}}{\rho_{\mathrm{h}}^{\max } \times \mathrm{SDI}_{i}{ }^{\max }}, $

(4)

其中：κ=0.5，μ=0.8，λ=0.5为以上3种因素对站域指数SI_i的贡献因子；ρ_h^max为静态设备随时间变化的流量比例的最大值；SDI_i^max为静态设备密度指数最大值。当站域指数SI_i(SI_i>0)大于阈值ω_SI(0 < ω_SI < 100%)时则认为该区域为站域，否则为非站域。

由于站域卫星存在比较大的拥塞可能性，而端到端流量无需经过地面站进行星地传输，因此在计算路径时尽可能少地使用站域卫星作为中间节点，减轻站域的流量负载。为了区分不同链路的状态并实现分流，本文根据站域指数SI_i分区计算卫星S_i与卫星S_j之间链路代价cost^ij，链路代价cost^ij由链路队列排队代价cost_que^ij及链路传播代价cost_prop^ij共同决定：

$ \operatorname{cost}^{i j}=\operatorname{cost}_{\text {que }}^{i j}+\operatorname{cost}_{\text {prop }}^{i j}, $

(5)

其中链路传播代价cost_prop^ij为链路传播时延T_prop^ij，即

$ \operatorname{cost}_{\text {prop }}^{i j}=T_{\text {prop }}^{i j}, $

(6)

$ T_{\text {prop }}^{i j}=\frac{d_{i j}}{c}, $

(7)

式中：d_ij为2颗卫星间距离，c为光速。

链路队列排队代价cost_que^ij主要由链路队列排队时延T_que^ij决定：

$ T_{\mathrm{que}}^{i j}=\frac{\mathrm{Q} \mathrm{OR}_{i j}}{\nu}, $

(8)

其中：QOR_ij为卫星S_i与临近卫星S_j链接的星间链路的队列占用率，ν为发送速率。当站域发生拥塞时，链路排队时延T_que^ij显著增大，进而引起链路队列排队代价cost_que^ij增加。为了研究站域卫星潜在拥塞可能性对链路排队时延T_que^ij的影响，引入站域潜在拥塞代价φ_i(φ_i≥0)来表示由站域卫星潜在拥塞可能性而带来的额外链路队列排队代价，对于非站域卫星不使用该参数，得到链路队列排队代价cost_que^ij计算如下

$ \mathrm{cost}_{\mathrm{que}}^{i j}=\left\{\begin{array}{c} T_{\mathrm{que}}^{i j}+\varphi_{i}, S_{i} \text { 为 } \mathrm{SAS}, \\ T_{\mathrm{que}}^{i j}, S_{i} \text { 为 } \mathrm{nSAS}. \end{array}\right. $

(9)

站域指数SI_i越高，则说明覆盖该区域的卫星越容易发生拥塞，故站域潜在拥塞代价φ_i计算如下

$ \varphi_{i}=S I_{i} \times T_{\text {que }}^{i j}, S_{i} \text { 为 } \mathrm{SAS}. $

(10)

综上所述，可以计算得到全网链路代价并作为多径计算过程中链路权重因子ψ_ij计算的基础。

2.3 多径搜索算法

TLR算法中使用最短路径算法进行多径计算得到一条最优路径及一条备选路径。本文提出的RMLBR算法在站域划分的基础上进行区域分流，尽可能选择链路代价cost^ij小的卫星作为中间节点，从而使较为空闲的卫星得到利用，缓解站域卫星的负担，实现负载均衡。此外，RMLBR算法以当前节点与目的节点之间距离d_id的倒数作为转向因子η_ij并引入节点可见性参数Γ_i，避免在多径搜索过程中出现绕行和环路。最后，为了满足在实时视频传输等多种时延敏感型场景中卫星数据的时效性，RMLBR算法以路径总时延T_path作为多径搜索的约束，进而保证数据传输的时效性。综上所述，RMLBR算法使用链路权重因子ψ_ij、转向因子η_ij及节点可见性参数Γ_i计算多径搜索中由节点S_i到节点S_j的转移概率p_ij，选择转移概率p_ij最大且满足约束的节点进行多径搜索进而得到最优路径及备选路径。具体过程详述如下：

首先，RMLBR算法使用链路权重因子ψ_ij表示链路代价cost^ij对路径计算的影响。链路代价cost^ij描绘网络中各条链路及各区域的状态，链路代价cost^ij越大时，则说明该链路不适合用于数据传输，进而算法选择经过该链路的概率越低。由此本文使用链路代价costij计算链路权重因子ψ_ij，即

$ \psi_{i j}=\omega / 1+\operatorname{cost}^{i j}, $

(11)

其中: ω为常量，本文中取ω=1.

其次，由于仅考虑链路代价而计算路径时容易舍近求远发生绕行，为减少冗余的中间节点并选择靠近目的节点的卫星作为中间跳，RMLBR算法以当前节点S_i与目的节点S_d之间的距离的倒数作为转向因子η_ij，即

$ \eta_{i j}=1 / d_{i d} $

(12)

此外，为避免在计算路径时出现环路，RMLBR算法设置了节点可见性参数Γ_i来标记该节点是否已被访问：

$ \varGamma_{i}=\left\{\begin{array}{l} 0, S_{i} \text { 已被访问, } \\ 1, S_{i} \text { 末被访问. } \end{array}\right. $

(13)

根据上述分析，RMLBR算法定义由节点S_i到临近节点S_j的转移概率p_ij为

$ p_{i j}=\left\{\begin{array}{c} \frac{\varGamma_{i} \times\left[\psi_{i j}\right]^{\alpha} \times\left[\eta_{i j}\right]^{\beta}}{\sum\nolimits_{S_{x} \in N(i)}\left[\psi_{i j}\right]^{\alpha} \times\left[\eta_{i x}\right]^{\beta}}, S_{j} \in N(i) ,\\ 0, S_{j} \notin N(i). \end{array}\right. $

(14)

其中: ψ_ij为链路权重因子；η_ij为转向因子；Γ_i为节点可见性参数；α，β为链路权重因子ψ_ij与转向因子η_ij的贡献系数；N(i)为当前节点S_i的临近节点集合。

在多径搜索策略开始时，将全网所有节点的可见性参数Γ_i均设置为1，表示所有节点均未被访问，并根据全网信息收集阶段得到的信息对链路权重因子ψ_ij及转向因子η_ij进行初始化。

最后，由于端到端时延是空间通信系统服务质量的重要影响因素^[15]，为保证时延敏感型场景下的服务质量及数据分组时效性，RMLBR算法引入路径时延约束。RMLBR算法将由源节点开始遍历后所经过的路径时延作为寻径约束，即已选择的节点组成的路径Pathm中各条链路的传播时延和排队时延之和T_path不能超过规定的路径时延的门限T_limit，其中路径时延T_path计算公式如下

$ T_{\text {path }}=\sum T_{\text {link }}^{i j}, $

(15)

$ T_{\text {link }}=T_{\text {que }}^{i j}+T_{\text {prop }}^{i j}. $

(16)

其中：T_link^ij为遍历过程中各条链路的时延。门限T_limit由实际场景中的时延要求决定，当门限T_limit取值过大时会失去约束作用，而取值过小易导致多径计算的失败。

综上所述，RMLBR算法由当前节点S_i选择临近节点S_j进行遍历需满足以下条件

$ \left\{\begin{array}{c} {\max}p_{i j}=\left\{\begin{array}{c} \frac{\varGamma_{i} \times\left[\psi_{i j}\right]^{\alpha} \times\left[\eta_{i j}\right]^{\beta}}{\sum\nolimits_{S_{x} \in N(i)}\left[\psi_{i j}\right]^{\alpha} \times\left[\eta_{i x}\right]^{\beta}}, S_{j} \in N(i), \\ 0, S_{j} \notin N(i). \end{array}\right. \\ \text { s. t. } T_{\text {path }} \leqslant T_{\text {limit }}. \end{array}\right. $

(17)

其中：α，β为链路权重因子与转向因子的贡献系数，N(i)为当前节点S_i的临近节点集合。

文献[16]提出，当多径之间存在公共节点时，其链路失效的可能性较大，因此为保证多径算法的性能，RMLBR算法中多径计算遵循最优路径及备选路径之间无公共节点的原则。在备选路径计算中屏蔽最优路径中选择的中间节点，使得构成最优路径的节点在备选路径计算中不可见，以此实现2条路径间无公有节点。RMLBR算法中的多径搜索策略具体步骤如下：

步骤1   对全网所有节点进行初始化：将卫星集合S中所有节点S_i可见性参数Γ_i均设置为1，并根据式(11)及式(12)对链路权重因子ψ_ij及转向因子η_ij进行初始化。

步骤2   搜索当前节点S_i的下一跳。若当前节点S_i为目的节点则算法结束；反之，根据式(14)计算转移概率p_ij并从节点S_i的临近节点集合N(i)中选择概率最大的节点S_j作为下一跳并转至步骤3。当节点S_j∈N(i)对应可见性参数Γ_j均为0时，说明N(i)中节点均已被访问，则退至节点S_i的前一跳节点S_pre重复步骤2进行搜索。

步骤3   对于已选择的下一跳节点S_j，计算从节点S_i到节点S_j的链路时延T_link^ij，若T_path+T_link^ij>T_limit，则节点S_j不满足约束，退回到节点S_i并转至步骤2重新进行搜索；若T_path+T_link^ij≤T_limit，则更新T_path=T_path+T_link^ij，Γ_j=0；并对节点S_j执行与节点S_i相同的搜索操作直至到达目的节点。

步骤4   在得到最优路径之后，重复步骤2~3计算备选路径。最终得到2条路径，即路径数目Num_path=2。

此外，RMLBR算法中多径搜索策略伪代码如表 1所示。

2.4 多径转发策略

在多径搜索结束后，全网各节点均将最优路径及备选路径写入路由表用于后续转发。由于在数据分组转发的过程中网络状态不断变化，因此需要在转发过程中根据当前的网络状态对路径做出调整，选择合适的下一跳节点进行分流，从而实现负载均衡。

RMLBR算法采用文献[11]中的“红绿灯”策略选择下一跳，根据卫星S_i各条链路的队列占用率QOR_ij及卫星S_i整体队列占用率TQOR_i设置卫星S_i各个方向上的红绿灯状态。红绿灯为“绿色”表示该方向未发生拥塞，“黄色”表示该方向即将拥塞，“红色”表示该方向已发生拥塞。当数据分组到达卫星S_i时，先从路由表中得到下一跳的候选，然后判断最优路径及备选路径中下一跳方向上的红绿灯状态选择合适的转发方式，规则如下：

1) 当最优路径下一跳方向上的红绿灯状态为“绿色”时，无论备选路径方向上是何种状态均选择最优路径下一跳进行转发。

2) 当最优路径下一跳方向上的红绿灯状态为“黄色”时，若备选路径下一跳方向上红绿灯状态为“绿色”或“黄色”，则进行分流，一半数据分组使用最优路径下一跳进行转发，另一半使用备选路径下一跳进行转发；若为“红色”，则选择最优路径下一跳进行转发。

3) 当最优路径下一跳方向上的红绿灯状态为“红色”时，若备选路径下一跳方向上红绿灯状态为“绿色”或“黄色”，则使用备选路径下一跳进行转发；若为“红色”，则令数据分组在缓存区等待至任一路径为非红色状态再进行转发。

3 算法仿真 3.1 仿真参数设置

为了研究本文提出的RMLBR算法在具有地面站的低轨卫星通信系统中的性能，利用OPNET modeler进行网络仿真，以中国科学院某在研类铱星极轨道试运行阶段的星座作为仿真背景，其中地面站位置根据该星座设计说明设置。在多种网络总流量输入下，从高流区域形态、网络丢包率、端到端时延及网络吞吐量4个方面与TLR算法和ELB算法进行对比，分析算法特性。

此外，本文不考虑信道误码对系统性能产生的影响，仅分析路由层算法性能。仿真过程中将上文所述流量模型作为参考输入进行研究，全网星间链路均配置200 Mbps容量，各个队列容量均为66 Mbit。并根据实际系统应用场景及服务质量需求将RMLBR算法中站域指数阈值ω_SI设置为0.6，路径时延门限T_limit设置为280 ms，ELB算法及TLR算法中的相关参数完全按照文献[10-11]设定。由于文中流量模型周期为24 h，故设置仿真时间为24 h并多次运行取均值进行研究。

3.2 仿真结果及分析

根据文献[4]，本文将卫星队列占用率TQOR_i>40%的区域视为高流区，图 5中阴影部分表示在网络总流量输入为148.3和189.5 Tbit/d时，3种算法的实时高流区形态，其中RMLBR算法的高流区范围最小，表明该算法有着较好的分流性能且缓解了拥塞。站域的划分预测了潜在的易拥塞区域，链路代价的分区计算能够引导路径搜索时尽可能少地使用站域卫星作为中间节点，降低了站域卫星的负载，使得较为空闲的卫星得到使用，从而实现负载均衡。

图 6展示多种不同的全网总流量输入下，3种算法平均端到端时延、丢包率及网络吞吐量情况。从图 6(a)可以看出，RMLBR算法的平均端到端时延最低。其中，在148.3 Tbit/d的输入下，ELB、TLR、RMLBR算法的端到端时延分别为118.8、110.9和107.3 ms，RMLBR算法相对降低了端到端时延，其原因在于该算法在多径计算时加入了时延的约束并使用转向因子以减少路径中出现绕行的情况。端到端时延的缩短使RMLBR算法在实时视频传输等多种时延敏感型场景下具有更高的实用性。

此外，由图 6(b)和6(c)可见，RMLBR算法在丢包率及吞吐量方面也具有优势，以148.3 Tbit/d输入为例，ELB、TLR、RMLBR算法的丢包率分别为6.68%、5.03%和3.38%，吞吐量分别为790.41、811.51和830.51 Gbit。RMLBR算法的性能优势主要得益于算法的设计充分考虑站域带来的潜在拥塞可能性计算端到端链路代价，并进行区域分流，从而缩小高流区范围，缓解了星地传输引起的空间段动态漏斗型拥塞。

综上所述，通过对ELB、TLR、RMLBR算法的高流区域形态、网络丢包率、端到端时延及网络吞吐量4个方面的仿真验证了RMLBR算法的有效性。

4 总结

本文针对低轨卫星通信系统非均匀分布地面信关站就近传输引起的空间段动态漏斗型拥塞问题进行研究，提出一种适用于对时延要求较高的场景下的路由算法——RMLBR。仿真结果表明，相对于经典的ELB算法和TLR算法，该算法能够缩小高流区域范围，缓解站域拥塞，并有效地降低端到端时延及网络丢包率，提高了网络吞吐量。