注水法在流量工程（Traffic Engineering）中的应用

本文链接：https://blog.youkuaiyun.com/lyhaoyuhangyuan/article/details/146148375

1. 核心思想

注水法的核心理念是将流量优先分配到资源更充裕（如带宽大、延迟低、拥塞程度轻）的路径或链路上，类似于注水时水自然流向低洼处。其目标是通过优化流量分配，最小化链路利用率、降低延迟或平衡负载，从而提升网络效率。

2. 数学模型

优化问题建模

典型问题可建模为凸优化问题，例如：

目标函数：最小化总延迟（基于M/M/1模型）
$\min \sum_{l} \frac{f_l}{c_l - f_l}$
其中， $f_{l}$ 为链路 $l$ 的流量， $c_{l}$ 为链路容量。
约束条件：
- 流量守恒： $\sum_{i} x_{i}^k = d^k$ （需求 $d^{k}$ 的流量需被满足）。
- 链路容量限制： $f_l = \sum_{k} x_{i}^k \leq c_l$ 。
- 非负性： $x_{i}^{k}\geqslant 0$ 。

拉格朗日对偶分解

通过引入对偶变量（链路价格 $\lambda _{l}$ ），构造拉格朗日函数：
$L = \sum_{l} \frac{f_l}{c_l - f_l} + \sum_{l} \lambda_l (f_l - c_l) +$ 流量守恒约束项.
对偶变量 $\lambda _{l}$ 反映链路的拥塞成本，优化过程需调整流量分配 $x_{i}^{k}$ 和对偶变量，直至达到均衡。

3. 算法步骤

初始化：设置链路价格 $\lambda _{l}^{(0)}$ 和初始流量分配 $f_{l}^{(0)}$ 。
流量分配：根据当前 $\lambda _{l}$ ，各节点选择路径使边际成本（延迟增量）等于路径价格总和。
- 路径价格： $\mu _{k}=\sum_{l\in k}\lambda _{l}$ 。
- 分配原则：将流量分配至 $\mu _{k}$ 最低的路径，直至各路径边际成本相等。
更新对偶变量：根据链路负载调整价格：
$\lambda _{l}^{(t+1)}=[\lambda _{l}^{(t)}+\alpha(f_{l}-c_{l}) ]^{+}$ ,
其中 $\alpha$ 为步长， $[.]^{+}$ 表示非负约束。
迭代收敛：重复步骤2-3，直至流量分配和对偶变量稳定。

4. 应用示例

场景设定

假设某网络流量需求为 $d=10 Gbps$ ，需通过两条路径传输：

路径1：容量 $c_{1}=10 Gbps$ ，当前负载 $f_{1}=7 Gbps$ 。
路径2：容量 $c_{2}=8 Gbps$ ，当前负载 $f_{2}=3 Gbps$ 。

目标是通过注水法重新分配流量，使得两条路径的链路价格（拥塞成本）相等，以最小化总延迟。

步骤详解

1. 初始状态

路径1：剩余容量 $c_{1}-f_{1}=3 Gbps$ ，链路价格：
$\lambda _{1}=\frac{1}{(c_{1}-f_{1})^{2}}=\frac{1}{3^{2}}=0.111.$
路径2：剩余容量 $c_{2}-f_{2}=5 Gbps$ ，链路价格：
$\lambda _{2}=\frac{1}{5^{2}}=0.040.$

此时路径2的链路价格更低，说明其拥塞成本更小，应优先分配流量至路径2。

2. 流量重新分配

根据注水法原则，流量应从高价格路径（路径1）转移到低价格路径（路径2），直到两路径价格相等。
设从路径1转移 $\Delta Gbps$ 到路径2，则新负载为：

路径1： $f_{1}^{'}=7-\Delta$
路径2： $f_{2}^{'}=3+\Delta$

需满足约束：

路径1负载不超过容量： $7-\Delta \leq 10$ （恒成立）。
路径2负载不超过容量： $3+\Delta \leq 8\Rightarrow \Delta \leq 5.$

3. 求解均衡条件

当两路径价格相等时，满足：

$\frac{1}{(10 - (7 - \Delta))^2} = \frac{1}{(8 - (3 + \Delta))^2}$

化简方程：

$\frac{1}{(3+\Delta )^{2}}=\frac{1}{(5-\Delta )^{2}}\Rightarrow 3 + \Delta = 5 - \Delta \Rightarrow \Delta = 1.$

4. 更新流量分配

转移 $\Delta =1Gbps$ 后：

路径1： $f_{1}^{'}=7-1=6Gbps$
剩余容量 $10-6=4Gbps$ ，新链路价格：
$\lambda_1' = \frac{1}{4^2} = 0.0625$
路径2： $f_{2}^{'}=3+1=4Gbps$
剩余容量 $8-4=4Gbps$ ，新链路价格：
$\lambda_2' = \frac{1}{4^2} = 0.0625$

此时两路径价格相等（ $\lambda _{1}^{'}=\lambda _{2}^{'}$ ），达到最优分配。

5. 结果验证

总延迟计算（基于M/M/1模型）：
- 路径1延迟： $\frac{6}{10 - 6} = 1.5$ 单位时间
- 路径2延迟： $\frac{4}{8 - 4} = 1.0$ 单位时间
- 总延迟： $1.5+1.0=2.5$ 单位时间
对比初始状态：
- 初始路径1延迟： $\frac{7}{3}\approx 2.333$
- 初始路径2延迟： $\frac{3}{5}= 0.6$
- 初始总延迟： $2.333+0.6\approx 2.933$