部分观测肾图模型中动力学速率常数在含噪时间-活性曲线下的实际可恢复性与数值可辨识性研究

《Molecular Imaging》：Practical Recoverability and Numerical Identifiability of Kinetic Rate Constants from Noisy Time-Activity Curves in a Patient-Derived Renogram Model

　　

摘要专业翻译 模型化时间-活性曲线（TACs）与观测数据之间的高度吻合，并不必然意味着潜在动力学参数能被准确恢复。研究人员利用一个具有潜在储库（latent reservoir）且仅观测左肾、右肾及膀胱TACs的部分观测四房室肾图模型，系统评估了上述区别。研究基于六份去标识化肾图数据提供患者来源的参考参数向量，开展添加噪声的合成恢复实验，并结合时间采样、优化器行为、初始化及数值可辨识性分析。相对于干净合成TACs的中位总归一化均方根误差（NRMSE）约为10?3量级，表明TAC重建效果优异。但参数恢复的均匀性明显较差：部分速率常数可被准确恢复，而另一些则出现数十至数百倍的乘性误差、反复触及边界，或涉及潜在储库的补偿效应。重复优化亦稳定收敛至错误的参数向量，却仍保持优异的TAC拟合质量。缩放残差雅可比矩阵（scaled residual Jacobian）数值上满秩，但严重病态（ill-conditioned），且若干参数的95%剖面似然区间（profile likelihood intervals）在优化边界处仍呈开放状态。上述发现表明，仅凭TAC拟合质量不足以评估参数可靠性。TAC重建、参数恢复、数值秩、边界行为、优化器稳定性及剖面似然提供了互补证据，在评估部分观测动力学模型时应综合解读。

论文解读文章

一、研究背景与问题

在动态核医学研究中，时间-活性曲线（time-activity curves, TACs）以紧凑形式描述示踪剂在体内的动态分布过程。房室模型（compartment model）利用这些曲线估计转移速率与清除速率，其参数比经验曲线参数更具生理可解释性。然而，参数估计面临多重固有困难：TACs数据有限、含噪声且常为稀疏采样，动力学参数可能跨越多个数量级，并通过强耦合通路共同影响可观测房室。

结构可辨识性（structural identifiability）与实用可恢复性（practical recoverability）之间存在本质区别。前者探讨在理想无噪声观测条件下模型参数是否可被唯一确定；后者则考察在有限采样、测量噪声、宽初始化范围、参数边界及数值优化等现实条件下，参数能否被稳定推断。已有研究通过解析和图论方法建立了若干房室模型的可辨识性结论，但这些结果并不自动适用于本研究的观测结构——储库为潜在状态，仅左肾、右肾和膀胱TACs被测量。

此前基于优化的四房室框架已证明生理速率常数可拟合患者TACs并重建观测曲线，但遗留的关键问题是：当数据生成参数向量已知时，拟合得到的参数向量能否被准确恢复。由于实测患者TACs无法区分参数误差与生理异质性、模型偏差或产生近似相同输出的替代参数化方案，必须开展已知真值的直接实验。

本研究针对四个相关问题展开：混合泊松-高斯TAC噪声如何传播至各速率常数；采样密度与采集时长如何改变参数特异性恢复及近边界行为；恢复结论对优化器工作流、重复起点及初始化宽度的敏感程度；数值轨迹秩、雅可比矩阵条件数与剖面似然如何与含噪条件下的经验恢复相关联。

二、研究设计与主要技术方法

研究采用正向-逆向（forward–inverse）设计。首先将四房室模型拟合至六份去标识化动态肾图数据，所得参数向量作为固定数据生成输入，用于生成干净合成TACs，再经扰动、重采样并重新拟合，以量化相对已知模拟输入的参数恢复精度。六份参考研究覆盖不同TAC形态：患者1–3为主参考配置，S1–S3为预设压力测试，分别侧重边界主导行为、大储库初始活性（A_R0）补偿及两者组合。

模型为线性时不变四房室系统，状态向量包含潜在储库、左肾、右肾和膀胱活性。观测矩阵仅保留左肾、右肾和膀胱三个房室。七个月速率参数进入系统矩阵，A_R0通过初始条件进入。所有速率在logarithmic参数空间优化以保证正值并改善跨数量级参数的条件数。速率公共边界为10^?4至10 min^?1。A_R0边界由处方活性信息经物理衰变校正后减去初始观测房室活性得到，确保模型初始化时总活性不超过处方预算。

基线拟合工作流从38种配置网格中选定，结合CMA-ES（协方差矩阵自适应进化策略，种群16，初始步长0.25，最多1600次迭代）与有界非线性最小二乘抛光（最多10,000次函数评估）。每份参考研究完成10次多起点拟合，保留目标函数值最低的解。

参数恢复实验采用混合泊松-高斯噪声模型：泊松分量代表计数依赖统计噪声，有效计数间隔q=25 s；加性高斯分量代表计数统计之外的残余变异，其标准差为对应泊松样本标准差的25%。每份参考向量生成500个独立噪声实现，共3000次合成恢复拟合。恢复拟合采用方差加权残差目标，并针对每个实现重新计算A_R0上界。

三、研究结果

参数恢复实验

在3000次合成恢复拟合中，TAC重建误差中位总NRMSE约为10^?3量级，表明模型可极好地重建含噪TACs。然而参数恢复远非均匀：部分速率常数恢复准确，另一些则出现数十至数百倍乘性误差、反复触及边界，或涉及潜在储库的补偿效应。重复优化稳定收敛至错误参数向量，同时保持优异TAC拟合，说明目标函数存在多个近等效的平坦区域。

噪声水平与时间采样扫描

噪声水平升高加剧参数恢复退化，但不同参数受影响程度差异显著。采样密度降低与采集时长缩短对参数特异性恢复产生不同影响：某些参数在稀疏采样下仍可恢复，而另一些则迅速失去可辨识性。边界接触行为随采样方案改变而改变，提示观测信息量对参数约束能力的直接作用。

优化器调优与对比

不同优化器工作流对参数恢复结论具有实质影响。部分优化器虽能获得优异TAC拟合，却系统性收敛至错误参数区域；另一些优化器则对初始化更敏感。固定噪声多起点实验表明，同一噪声实现下不同起点可收敛至不同参数向量，而TAC拟合质量几乎相同，进一步证实目标函数的多模态性。

初始化宽度与数值灵敏度

初始化宽度扩大显著恶化参数恢复，尤其对弱可辨识参数。缩放残差雅可比矩阵数值上满秩，但条件数严重恶化，表明参数方向之间存在近线性依赖。Fisher信息分析与此一致，显示若干参数方向的信息量极低。

剖面似然分析

多个参数的95%剖面似然区间在优化边界处仍呈开放状态，意味着这些参数在给定数据下无法被有限界定。边界开放与经验恢复误差模式高度一致，确认了部分参数仅受边界而非数据约束。

四、讨论与结论

研究综合TAC重建、参数恢复、数值秩、边界行为、优化器稳定性与剖面似然六类证据，系统刻画了部分观测潜状态模型中TAC一致性与参数恢复之间的差距。核心结论为：TAC拟合质量本身不足以评估参数可靠性。即使模型可完美重建观测曲线，部分速率常数仍可能无法被数据唯一确定，且优化器可能稳定收敛至错误参数向量。

该发现对动态核医学定量分析具有重要方法学意义：在报告房室模型参数时，应同时呈现参数恢复诊断指标，而非仅依赖拟合优度。剖面似然开放区间与雅可比矩阵病态性可作为弱可辨识性的直接警示信号。研究同时提示，采样方案设计与优化器选择对参数可恢复性具有实质性影响，需在协议设计阶段予以考虑。

研究人员将合成真值定义为固定参数向量生成的模拟数据集，不暗示对相应临床肾图的独立验证性生理表征。研究结论局限于部分观测四房室模型框架，但其方法论框架——结合已知真值合成实验、多维度诊断指标与剖面似然分析——可推广至其他部分观测动力学系统的参数可靠性评估。