恢复扩散策略中的隐藏奖励 Recovering Hidden Reward in Diffusion-Based Policies
用标量能量函数统一扩散模仿学习与逆强化学习,无需对抗训练即可恢复奖励信号
前置知识
去噪分数匹配 (Denoising Score Matching)
一种估计数据分布对数密度梯度(即分数函数 $\nabla_x \log p(x)$)的训练目标。给定干净样本 $x_0$ 和噪声扰动核 $q_\sigma(\tilde{x}|x_0)=\mathcal{N}(\tilde{x};x_0,\sigma^2 I)$,其目标等价于让网络 $S_\theta$ 预测噪声方向 $-\epsilon/\sigma$,从而避开了直接计算难处理的配分函数。
本文的能量网络 $E_\phi$ 通过自动微分得到分数 $S_\phi = -\nabla_a E_\phi$,其训练过程正是用去噪分数匹配完成的;理解这一点是看懂 ENERGYFLOW 训练目标(公式 15)的前提。
扩散策略 (Diffusion Policy)
一种将模仿学习策略参数化为条件去噪扩散模型的范式(Chi et al., 2023)。它学习一个噪声条件分数网络 $S_\theta(a_t, s, t)$,推理时从 $\mathcal{N}(0,\sigma^2(T)I)$ 出发,用概率流 ODE 迭代去噪得到专家动作,能很好地建模多模态动作分布。
ENERGYFLOW 直接继承了 Diffusion Policy 的 1D U-Net 主干网络与 ODE 采样流程,但把'直接回归噪声'改为'输出标量能量再求梯度';理解原版框架的边界条件能帮助读者抓住这一关键改动。
最大熵逆强化学习 (Max-Entropy IRL)
由 Ziebart et al. (2008) 提出的 IRL 范式,假设专家策略满足软最优性 $\pi_E(a|s)=\frac{1}{Z(s)}\exp(Q^*(s,a)/\alpha)$,其中 $Q^*$ 是软动作价值函数。该假设让奖励恢复问题可以用分布匹配的方式求解,但配分函数 $Z(s)$ 通常难以计算。
本文 Theorem 3.3 正是建立在该假设之上,证明在最大熵最优性条件下 $\nabla_a Q^*(s,a) = \alpha \cdot S^*(a,s)$,即分数函数就是软 Q 函数的梯度,这是全文核心等式的源头。
保守向量场 (Conservative Vector Field)
若一个向量场 $V:\mathbb{R}^d \to \mathbb{R}^d$ 满足 $V = \nabla \Psi$ 即它是某标量势能 $\Psi$ 的梯度,则称其为保守(或可积)的;等价条件是其 Jacobian 对称 $\nabla \times V = 0$。保守场保证路径无关与偏好传递性。
论文强制 $S_\phi = -\nabla_a E_\phi$ 使分数场天然保守。Theorem 3.6 证明这一约束将假设类的 Rademacher 复杂度从 $\mathcal{O}(\Lambda B\sqrt{d/n})$ 收紧到 $\mathcal{O}(\Lambda L/\sqrt{n})$,从而改善泛化。
研究动机
扩散策略(Diffusion Policy)虽然能很好地拟合多模态专家动作分布,但其训练目标本质上还是行为克隆——只模仿"做什么",没有显式建模"为什么"。当测试时遇到与演示分布有偏差的状态,仅靠似然匹配往往不能给出可靠的动作选择信号,从而限制了鲁棒性和外推能力。RoboMimic Transport、ToolHang 等接触密集任务的真实成功率也常常被这一瓶颈拖累。另一方面,虽然用强化学习给扩散策略加 reward 是常见补救思路,但现实机器人中手工设计 reward 既繁琐又易错;现有的逆强化学习方法如 GAIL/AIRL 走对抗路线存在训练不稳定和模式坍缩问题,而基于能量模型的 EBIL、NEAR 又需要在高维动作空间里做昂贵的 MCMC 采样估计配分函数,难以扩展。一个关键观察是:扩散策略里的去噪分数网络其实已经在学 $\nabla_a \log p(a|s)$,而最大熵框架下 $\nabla_a \log \pi_E(a|s)$ 与软 Q 函数梯度是线性关系,因此分数场里其实"藏着"专家奖励,但现有方法没有把这层联系利用起来。
本文的目标是本文提出一个名为 ENERGYFLOW 的统一框架,目标是同时解决三件事:(1) 保持甚至超越现有扩散策略的模仿学习性能;(2) 在不需要对抗训练或 MCMC 的前提下,自动从专家演示中恢复出可解释的奖励信号;(3) 借助该奖励信号做下游强化学习,使策略能在环境交互中持续改进。作者明确希望把"做模仿"和"恢复奖励"两件事统一在一个去噪目标里完成,让结构约束同时充当奖励提取的合法性条件和泛化的归纳偏置。
与已有工作不同的是,现有工作大致有三类切入角度,但都不完整:第一类是 Diffusion Policy 系,仅做行为克隆,没有奖励通道;第二类是基于对抗或 MLE 的 IRL(GAIL/AIRL、IQ-Learn、EBIL/NEAR),要么训练不稳定,要么需要昂贵采样,且大多与扩散模型解耦;第三类是 Implicit BC、EBT-Policy 等端到端能量策略,把能量当作决策分数而非可识别奖励。ENERGYFLOW 的独特切入点是:通过直接对标量能量函数 $E_\phi(a,s,t)$ 而不是向量噪声做去噪分数匹配,利用 $\nabla_a E_\phi$ 的自动微分自动满足保守场假设,从理论上证明这等价于在最大熵假设下恢复软 Q 函数梯度,并由此得到一个零额外成本、低方差的标量 reward。
核心方法
ENERGYFLOW 的核心直觉非常简洁:与其让网络直接输出噪声 $\epsilon_\theta$(也就是分数场),不如让它先输出一个标量"能量"$E_\phi(a,s,t)$,然后用 PyTorch 的自动微分对动作求梯度 $\nabla_a E_\phi$,把梯度取负号就得到了 score。这种参数化在数学上强制了 $\nabla \times (\nabla_a E_\phi) \equiv 0$,即学到的向量场天然是保守的;与此同时,这块"保守约束"在最大熵最优性条件下让 $E_\phi$ 恰好恢复软 Q 函数(差一个状态相关常数)。训练仍沿用方差爆炸 (VE) 形式的去噪分数匹配,损失为 $\mathcal{L}=\mathbb{E}_{t,a_0,\epsilon}[\sigma^2(t)\|S_\phi(a_t,s,t)+\epsilon/\sigma(t)\|^2]$,其中 $S_\phi = -\nabla_{a_t}E_\phi$;噪声调度 $\sigma(t)=\sigma_{\min}^{1-t/T}\sigma_{\max}^{t/T}$,$\sigma_{\min}=0.01$,$\sigma_{\max}=10.0$,$T=1.0$。推理时和扩散策略一样用 Euler 方法沿概率流 ODE 从 $t=T$ 积分到 $t\approx 0$,但每一步的速度场换成了 $\frac{1}{2}\frac{d\sigma^2(t)}{dt}\nabla_a E_\phi$。下游使用时,作者提出"中心化塑形 (centered shaping)":$\tilde{r}_\phi(a,s) = -\big(E_\phi(a,s,\epsilon) - \mathbb{E}_{a'\sim\mathcal{N}(0,I)}[E_\phi(a',s,\gamma)]\big)$,用 $M=16$ 个参考样本估计基线,把 state-dependent 偏移减掉,只保留相对偏好。
本工作和已有方法的本质区别有两点。第一个区别是把"显式能量 + 自动微分梯度"作为结构约束,而不是把能量当作事后优化目标:IBC、EBT-Policy 把能量当作非凸黑盒靠 Langevin 动力学或 DFO 搜索最小值,没有理论保证恢复出 reward;ENERGYFLOW 通过让网络只输出标量、由 $\nabla_a$ 自动产生向量场,从架构层面就保证了保守性,无需额外正则。第二个区别是理论层面的"等价性证明":Theorem 3.3 严格说明在最大熵假设下 $\nabla_a Q^* = \alpha S^*$,由此学到的 $E_\phi$ 恢复出软 Q 函数(差一个状态相关常数);Theorem 3.6 和 Lemma 3.8 进一步把保守性约束转化为 Rademacher 复杂度从 $\mathcal{O}(\sqrt{d})$ 降到 $\mathcal{O}(L)$ 的泛化优势,并给出 OOD 风险界 $R_{D_T}(\hat{f}_{\text{cons}}) \leq \hat{R}_S(\hat{f}_{\text{cons}}) + \frac{1}{2}d_{\mathcal{H}\Delta\mathcal{H}}(D_S,D_T) + \mathcal{O}(\frac{M\Lambda L}{\sqrt{n}})$,而无约束场随 $d$ 发散。这意味着 ENERGYFLOW 不是事后给 diffusion policy 加一个 reward 通道,而是把"结构合法性"直接融进了训练目标本身。
方法步骤详情
完整流程包括三个阶段。第一阶段是网络架构与改造:在 Diffusion Policy 的 1D Conditional U-Net 基础上,把状态用 2 层 MLP(hidden=128, Mish)编码、时间用 sinusoidal embedding 后线性投影、两者通过 FiLM 注入每个卷积块;与原版最大的不同是输出头——用 GlobalAveragePooling1D 把时序特征聚合,再过 3 层 MLP(256→128→1)输出一个标量 $E_\phi$,并把所有 ReLU 换成 Mish 以保证 $C^2$ 可微(因为训练需要对 score 再求一次梯度)。第二阶段是去噪分数匹配训练:对每个 batch 采样 $(s,a_0)\sim \mathcal{D}$、$t\sim U[0,T]$、$\epsilon\sim\mathcal{N}(0,I)$,构造 $a_t = a_0 + \sigma(t)\epsilon$,调用 `torch.autograd.grad(..., create_graph=True)` 求 $S_\phi(a_t,s,t) = -\nabla_{a_t} E_\phi$,然后最小化 $\mathcal{L}=\sigma^2(t)\|S_\phi(a_t,s,t) + \epsilon/\sigma(t)\|^2$;优化器使用 AdamW(lr=$10^{-4}$,wd=$10^{-6}$,batch=256,Cosine Decay,warmup=500,梯度裁剪 norm=1.0),并对输出头线性层做谱归一化以控制 Lipschitz。第三阶段是推理与 reward 提取:动作生成时(Algorithm 2)以 $a_T \sim \mathcal{N}(0,\sigma^2(T)I)$ 为起点,循环 $K=20$ 步执行 Euler 更新 $a_{k+1} = a_k - \Delta t \cdot g_k$,其中 $g_k = \frac{1}{2}\frac{d\sigma^2(t_k)}{dt}\nabla_a E_\phi(a_k,s,t_k)$,最终返回 $a_K$ 与对应能量 $E_\phi(a_K,s,\gamma)$;奖励信号则按公式 (16) 做中心化塑形,用 $M=16$ 个 $a'\sim\mathcal{N}(0,I)$ 估计基线项 $\mathbb{E}[E_\phi(a',s,\gamma)]$,再用 SAC 在 200k 步环境交互中训练下游策略。
技术新颖性
技术新颖性可以从三个层次看。架构层,把"标量能量 + 自动微分梯度"引入扩散策略,配套的改造——全局时序池化 + 标量 MLP 输出头 + Mish + `create_graph=True` 的双反传——是为了同时满足 $C^2$ 可微与保守性约束,这些都不是简单把 Diffusion Policy 改个 head 就能做到的。理论层,论文给出了完整的"score-reward 等价"证明链:Assumption 3.1(最大熵最优) → Theorem 3.3($\nabla_a Q^* = \alpha S^*$) → Corollary 3.4(恢复 soft advantage) → Proposition 3.9(同状态内 ranking 精确,跨状态有常数模糊) → Theorem 3.6(保守类 Rademacher 复杂度 $\leq \Lambda L/\sqrt{n}$ 而非约束类 $\leq \Lambda B\sqrt{d/n}$) → Lemma 3.8(OOD 风险界随 $d$ 控制而非发散) → Theorem 3.11(score 估计误差 $\epsilon$ 线性传递到偏好误差 $\epsilon\|a-a'\|_2$),这套证明在 diffusion policy 领域比较少见。系统层,"同一目标同时支撑模仿学习 + 奖励提取 + 下游 RL"是端到端流水线,相比 EBIL/NEAR 那种"先学能量、再 RL"的两阶段方案,少了模式坍缩和配分函数估计的麻烦。
实验结果
论文从 6 个研究问题展开实验,每个都有具体数字支撑。RQ1(模仿性能):RoboMimic 五任务平均 ENERGYFLOW 取得 93.8% 成功率,超过 Diffusion Policy 的 91.2% 和 Flow Policy 的 89.6%;其中最难的 ToolHang 任务提升最大,从 77.2% 提升到 84.2%,Lift 和 Can 都达到 100%。Meta-World 五任务平均 92.5%,超过 Diffusion Policy 的 90.7%;现有 energy-based 方法如 EBT-Policy 平均只有 77.8%,Implicit BC 在 RoboMimic Can 上仅 30.8%,更复杂任务直接 0.0%。RQ2(真实机器人):用 AGIBOT G1(7-DoF + 平行夹爪,单目 RGB)在 Bottle(抓瓶子放进纸箱)和 Drawer(拉开抽屉)两个接触密集任务上各采集 20 条专家演示,3 个初始位置扰动 × 20 rollout 都取得 100% 成功率,且轨迹比 Diffusion Policy 更平滑、接触点附近犹豫更少。RQ3(reward 质量):在 RoboMimic Square 和 Transport 上用 SAC 训练 200k 步对比四种 reward 信号,centered energy 取得接近 oracle dense reward 的成功率,而 raw energy 会陷入"停留在常见状态"的早期 plateau,sparse reward 早期几乎没信号;centered + sparse 组合效果最佳,说明 energy 提供逐步引导、sparse 兜底任务完成。RQ4(OOD 鲁棒性):在 RoboMimic 五个任务上施加 S/M/L 三级初始位置扰动,所有方法在 in-distribution 都接近饱和,但随着扰动加大,ENERGYFLOW 的优势扩大,在 L 级别扰动下保持显著高于 Diffusion Policy / Flow Policy 的成功率,验证了 Lemma 3.8 的预测。RQ5(reward 提取时间 $\gamma$ 敏感性):$\gamma\in[10^{-4},10^{-2}]$ 三档都在 88% 以上,$\gamma=10^{-3}$ 时最高 95.3%,$\gamma\geq 0.1$ 时因噪声分布偏离数据分布而掉到 72.6%,说明 score 近似只在低噪声段有效。RQ6(推理效率):在 NVIDIA A100 上 10Hz 控制,ENERGYFLOW (K=20) 延迟 11.4ms、成功率 95.3%,与 Flow Policy 的 8.2ms / 91.8% 相当,但比 Implicit BC (50 Langevin) 52.4ms / 10.2% 快得多;DDPM 100 步则要 32.4ms 且无能量输出。
查看结构化数据
| 任务 | 指标 | 本文 | 基线 | 提升 |
|---|---|---|---|---|
| RoboMimic (Lift/Can/Square/Transport/ToolHang 平均) | Success rate (%) | 93.8 | 91.2 (Diffusion Policy) | +2.6 |
| RoboMimic ToolHang | Success rate (%) | 84.2 | 77.2 (Diffusion Policy) | +7.0 |
| Meta-World (Button/Drawer/Assembly/Bin/Hammer 平均) | Success rate (%) | 92.5 | 90.7 (Diffusion Policy) | +1.8 |
| RoboMimic Can | Success rate (%) | 100.0 | 99.2 (Diffusion Policy) | +0.8 |
| Real Robot Bottle & Drawer | Success rate (%) over 60 rollouts | 100 | — (Diffusion Policy 仅定性对比) | 持平或更平滑 |
| RoboMimic Square (10Hz control) | Latency (ms) / Success (%) | 11.4 ms / 95.3% (K=20) | 9.1 ms / 90.4% (Diffusion Policy 20 DDIM); 52.4 ms / 10.2% (Implicit BC 50 Langevin) | 兼顾速度与可暴露标量能量 |
| Implicit BC 对比 (RoboMimic 平均) | Success rate (%) | 93.8 | 22.4 (Implicit BC) | +71.4 |
局限与改进
作者在 Remark 3.10 与 Theorem 3.11 附近诚实指出了几个边界:(1) 学到的 reward $\hat r = -\alpha E_\phi$ 与真实 soft Q 只差一个状态相关常数 $c(s)$,Proposition 3.9 证明同状态内 ranking 完全精确,但跨状态比较存在 $c(s)-c(s')$ 的模糊;在 sequential MDP 中这不一定满足 potential-based reward shaping 形式,可能改变最优策略,因此更适合作为 shaping / within-state 排序信号,而非直接 ground-truth reward。(2) Theorem 3.11 要求学到的分数满足全局 $\epsilon$ 误差界,最大单步 reward 估计误差为 $\alpha \cdot \epsilon \cdot \mathrm{diam}(\mathcal{A})$,在动作空间非常大或 score 拟合不精确时会放大。(3) 论文假设 Lipschitz 连续性靠权重衰减与谱归一化保证(Remark 3.7),但实际训练仍可能因双反传 (`create_graph=True`) 出现显存与稳定性开销。从我自己的观察看,还有一些隐含限制:仿真任务的状态是低维 proprioceptive(关节角/速度/物体位姿),没有覆盖 3D Diffusion Policy 那类视觉输入场景;real robot 仅 2 个任务、20 条演示,统计显著性有限;仿真平均成功率虽最高,但单任务如 Meta-World Drawer 仍输给 Diffusion Policy 0.6 个百分点(94.2 vs 93.6),优势主要集中在较难任务(ToolHang、Transport),说明"能量约束"对高维复杂分布建模有正面影响,但对简单任务不一定稳定占优。
独立分析的弱点
独立看几个具体弱点和对应改进方向。第一,$C^2$ 可微要求带来工程负担:必须把 ReLU 全部换成 Mish、GlobalAveragePooling 后才能用标量 MLP 头,导致时序局部信息被强制压缩成单值,理论上可能丢失精细时间结构;改进思路是探索局部能量参数化(如对每个时间步输出一个能量再相加)以兼顾局部细节与全局保守性。第二,state-dependent 常数 $c(s)$ 让 energy 不能直接当作 absolute reward,作者用 centered shaping 减去 $\mathbb{E}_{a'\sim\mathcal{N}(0,I)}[E_\phi(a',s,\gamma)]$ 来缓解,但该基线分布未必匹配真实 in-distribution 分布,下游 SAC 仍依赖 sparse 信号兜底;改进方向是引入 state-conditioned baseline 网络或用 demonstration state 分布做 importance sampling。第三,Theorem 3.6 的复杂度界要求表示 $\phi$ 满足 $L\ll B\sqrt{d}$,在欠平滑或高秩 Jacobian 的网络里保守约束的优势会被削弱;改进方向是在损失里加一个轻量的 symmetry penalty $\|\nabla_a\times S_\phi\|^2$(注:解析上恒为 0,可改为基于差分的有限差分近似),从而在无法严格 $C^2$ 时仍维持 curl-free 约束。第四,训练需 `create_graph=True` 的二阶反传,显存与算力开销大约是标准 Diffusion Policy 的 1.5–2 倍,论文没给出具体开销数字;改进方向是探索对偶数法或 Neumann 级数近似 Hessian。第五,OOD 实验只在 RoboMimic 的初始位置扰动上做了 S/M/L 三档,未覆盖视觉外观、动力学参数变化等更现实的分布偏移。
未来方向
作者在结论部分隐含了几条延伸方向。基于本文成果可以继续展开的研究包括:(1) 把能量恢复与离线强化学习(CQL、IQL)结合,利用 $E_\phi$ 作为保守 Q 函数初值,避免从头训练 critic;(2) 把 Theorem 3.3 的 score-reward 等价推广到 hierarchical / skill-conditioned diffusion policies,让组合策略也能恢复出组合 reward;(3) 把 centered shaping 与 uncertainty estimation 结合,给出能量估计的置信区间,做 safe exploration;(4) 探索把保守性约束与 Flow Matching 框架结合,因为 Flow Policy 已经在 inference 速度上优于 DDPM,理论上也能维持等价保守性而进一步提速;(5) 把 ENERGYFLOW 扩展到语言条件、多模态视觉输入的 setting,作者提到 3D Diffusion Policy 和 DexVLA 等方向都还没有显式 reward 通道,是自然的延伸;(6) 用 $E_\phi$ 做 reward shaping + curriculum learning,根据当前 state 的能量分布自动调整任务难度。
复现评估
复现评估总体良好但有一些门槛。优点是作者明确说代码开放在 https://github.com/sotaagi/EnergyFlow,Appendix C 给出了完整超参数(AdamW, lr=$10^{-4}$, wd=$10^{-6}$, batch=256, U-Net 通道 $[64,128,256]$, Mish, GroupNorm 8, Kernel 5, 谱归一化, Cosine Decay 500 warmup, 梯度裁剪 norm=1.0, ODE 步数 $K=20$, $T_p=16$, $T_o=2$),并公开了噪声调度 $\sigma(t)=\sigma_{\min}^{1-t/T}\sigma_{\max}^{t/T}$、centered shaping 的 $M=16$ 样本、真实机器人的 ResNet-18 ImageNet 预训练 backbone 等关键细节;附录 D 详述了 10 个仿真任务的内容。门槛主要是:(1) 训练需要二阶反传,单卡显存压力较大,作者虽提到 NVIDIA A100 但未报告 GPU 时长,按经验估计 RoboMimic 单任务训练约需 12–24 小时;(2) 真实机器人实验依赖 AGIBOT G1 + 10 Hz 控制 + ResNet-18,普通实验室复现成本较高;(3) 理论证明(Appendix A.1–A.4)虽然详尽,但实现 Theorem 3.6 中的 Jacobian Frobenius 范数界需要手工监控特征表示的 $L$ 与 $B$,对工程经验有要求;(4) 数据集方面,RoboMimic / Meta-World 都是公开 benchmark,但真实机器人的 20 条 Bottle / Drawer 演示没有公开,复现真实实验需要自行采集。总体判断:仿真部分可复现性 4/5(依赖算力),真实机器人部分 3/5。
论文图表
罗列架构(1D Conditional U-Net, Global Average Pooling, Mish)、训练(AdamW, lr=$10^{-4}$, wd=$10^{-6}$, batch=256, Cosine Decay 500 warmup, 梯度裁剪 1.0, 几何噪声调度)和推理(Euler, K=20, Tp=16, To=2)的全部超参。
复现必查的清单,Appendix C 的文字描述对应的总览表。